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This paper deals with the computation of sectional curvature for the manifolds of A'^ land- 
marks (or feature points) in D dimensions, endowed with the Riemannian metric induced by the 
group action of diffeomorphisms. The inverse of the metric tensor for these manifolds (i.e. the 
cometric), when written in coordinates, is such that each of its elements depends on at most 2D 
of the ND coordinates. This makes the matrices of partial derivatives of the cometric very 
sparse in nature, thus suggesting solving the highly non-trivial problem of developing a for- 
mula that expresses sectional curvature in terms of the cometric and its first and second partial 
derivatives (we call this Mario's formula). We apply such formula to the manifolds of landmarks 
and in particular we fully explore the case of geodesies on which only two points have non-zero 
momenta and compute the sectional curvatures of 2-planes spanned by the tangents to such 
geodesies. The latter example gives insight to the geometry of the full manifolds of landmarks. 

1 Introduction 

In the past few years there has been a growing interest, in diverse scientific communities, in modeling 
shape spaces as Riemannian manifolds. The study of shapes and their similarities is in fact central 
in computer vision and related fields (e.g. for object recognition, target detection and tracking, 
classification of biometric data, and automated medical diagnostics), in that it allows one to recognize 
and classify objects from their representation. In particular, a distance function between shapes 
should express the meaning of similarity between them for the application that one has in mind. One 
of the most mathematically sound and tractable methods for defining a distance on a manifold is to 
measure infinitesimal distance by a Riemannian structure and global distance by the corresponding 
lengths of geodesies. 

Among the several ways of endowing a shape manifold with a Riemannian structure (see, for 
example, [17, 18, 20, 25, 28, 30]), one of the most natural is inducing it through the action of 
the infinite-dimensional Lie group of diffeomorphisms of the manifold ambient to the shapes being 
studied. You start by putting a right-invariant metric on this diffeomorphism group, as described 
in [27]. Then fixing a base point on the shape manifold, one gets a surjective map from the group of 
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difFeomorphisms to the shape manifold. The right-invariance of the metric "upstairs" imphes that 
we get a quotient metric on the shape manifold for which this map is a submersion (see below). 
This approach can be used to define a metric on very many shape spaces, such as the manifolds 
of curves [12, 26], surfaces [33], scalar images [4], vector fields [6], diffusion tensor images [5], mea- 
sures [11, 13], and labeled landmarks (or "feature points") [14, 15]. The actual geometry of these 
Riemannian manifolds has remained almost completely unknown until very recently, when certain 
fundamental questions about their curvature have started being addressed [25, 26, 32]. 

Among all shape manifolds, the simplest case of the manifold of landmarks in Euclidean space 
plays a central role. This is defined as 

£^(M^) := {(Pi, . . . , P^) I P» e M^, a = 1, . . . , iv}. 

(typically we consider landmarks P", a = 1,...,N that do not coincide pairwise). It is finite- 
dimensional, albeit with high dimension n = ND, where N is the number of landmarks and D is 
the dimension of the ambient space in which they live (e.g. D = 2 for the plane). Therefore its 
metric tensor may be written, in any set of coordinates, as a finite-dimensional matrix. This space 
is important in the study of all other shape manifolds because of a simple property of submersions: 
for any submersivc map f : X ^ Y, all geodesies on Y lift to geodesies on X and give you, in fact, 
all geodesies on X which at one and hence all points are perpendicular to the fiber of / (so called 
"horizontal" geodesies). This means that geodesies on the space of landmarks lift to geodesies on the 
difi^eomorphism group and then project down to geodesies on all other shape manifolds associated to 
the same underlying ambient space M^. Thus geodesies of curves, surfaces, etc. in can be derived 
from geodesies of landmark points. Technically, these are the geodesies on these shape manifolds 
whose momentum has finite support. This efficient way of constructing geodesies on many shape 
manifolds has been exploited in much recent work, e.g. [2, 8, 29]. 

What sort of metrics arise from submersions? Mathematically, the key point is that the inverse 
of the metric tensor, the inner product on the cotangent space hence called the co-metric, behaves 
simply in a submersion. Namely, for a submersion / : X — >■ y, the co- metric on Y is simply the 
restriction of the co-metric on X to the pull-back 1-forms. Therefore, for the space of landmarks the 
cometric has a simple structure. In our case, we will see that each of its elements depends only on 
at most 2D of the ND coordinates. Hence the matrices obtained by taking first and second partial 
derivatives of the cometric have a very sparse structmc that is, most of their entries are zero. 
This suggests that for the purpose of calculating curvature (rather than following the "classical" 
path of computing first and second partial derivatives of the metric tensor itself, the Christoffcl 
symbols, et cetera) it would be convenient to write sectional curvature in terms of the inverse of 
the metric tensor and its derivatives. We have solved the highly non-trivial problem of developing 
a formula (that we call "Mario's formula") precisely for this purpose: for a given pair of cotangent 
vectors this formula expresses the corresponding sectional curvature as a function of the cometric 
and its first and second partial derivatives except for one term which requires the metric (but not 
its derivatives). This formula is closely connected to O'Neill's formula which, for any submersion as 
above, connects the curvatures of X and Y. Subtracting Mario's formula on X and Y gives O'Neill's 
as a corollary. 

This paper deals with the problem of computing geodesies and sectional curvature for landmark 
spaces, and is based on results from the thesis of the first author [23]. The paper is organized as 
follows. We first give a few more details about the manifold of landmarks, and describe the metric 
induced by the action of the Lie group of diffeomorphisms. We then give a proof for the general 
formula expressing sectional curvature in terms of the cometric. This formula is used in the following 
section to compute the sectional curvature for the manifold of labeled landmarks. In the last section, 
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wc analyze the case of geodesies on whieh only two points have non-zero momenta and the sectional 
curvatures of 2-planes made up of the tangents to such geodesies. In this case, both the geodesies 
and the curvature are much simpler and give insight into the geometry of the full landmark space. 



2 Riemannian Manifolds of Landmarks 

In this section we briefly summarize how the shape space of landmarks can be given the structure of 
a Riemannian manifold. We refer the reader to [27, 31] for the general framework on how to endow 
generic shape manifolds with a Riemannian metric via the action of Lie groups of diffeomorphisms. 



2.1 Mathematical preliminciries 

We will first define a distance function d : C^{R^) x £^(M^) R+ on landmark space which will 
then turn out to be the geodesic distance with respect to a Riemannian metric. Let Q be the set of 
differentiable landmark paths, that is: 



Q:={q={q\...,q^) : [0, 1] ^ £^(]R^) | € ([0, 1] , M^) 



L 
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Following [31, Chapters 9, 12, 13], a Hilbert space (V, ( , )y) of vector fields on Euclidean space 
(which we consider as functions — )■ M.^) is said to be admissible if (i) V is continuously embedded 
in the space of C^-mappings on R^ which are bounded together with their derivatives, (ii) V 

is large enough: For any positive integer M, if xi, . . . ,Xm € K'^ and ai, . . . , au & are such 
that, for all u eV, Y.^=i ("a: '^{^a))^o = 0, then ai = . . . = Um = 0. 

The space {V, { , )v) admits a reproducing kernel: that is, for each a,x E R^ there exists 
€ V with {KS,f}v = {a,.f{x)}^n for all / G V. Further, {K^,K^)v = {P,K^{v))^o = 
{a, K^{x))]g,n which is a bilinear form in (a, (3) € (M^)-^, thus given hy a. D x D matrix K{x, y); the 
symmetry of the inner product implies that K{y,x) = K{x,y)'^ (where indicates the transpose). 
In this paper we shall assume that K{x,y) is a multiple of the identity and is translation invariant: 
we then write K{x,y) simply as K{x — y)iD (where Id is the D x D identity matrix); the scalar 
reproducing kernel K : — )• R must be symmetric, and positive defijiite (see [31, §9.1] for details). 

There are other very natural admissible norms on vector fields v whose kernels are not multiples 
of the identity, e.g. one can add a multiple of div(i')^ to any norm and then K will intertwine 
different components of v. The most natural examples of the norms we will consider are given by 
inner products 

{u,v)v = {u,v)l := I {Lu{x),v{x)) o dx, (1) 

where L is a self-adjoint elliptic scalar differential operator of order greater than D + 2 with constant 
coefficients which is applied separately to each of the scalar components of the vector field u = 
{v^ , . . . ,u^). By the Sobolev embedding theorem then V consists of C^-functions on R^ which 
are bounded together with their derivatives, li K \s a. scalar fundamental solution (or Green's 
function [9]) so that L{K){x) = 5{x), then the reproducing kernel is given by K" = K{ —x) a. A 
possible choice of the operator is i = (1 — A^A)*^ (where A G M is a scaling factor, k G N and A is 
the Laplacian operator), with fc > y + 1, in which case (1) becomes the Sobolev norm: 



EE 

i=l m=0 



A' 



E 

|a|=m 



dx, 



(2) 
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When L = (1 — A^A)^ the scalar kernel K has the form K{x — y) = 7(||^ — 2/||mo)) with: 

1 k— — 

7(e) = , D . o (4) ' ^fc ^?>o, (3) 

where iiTi, (with v = k — ^) is a modified Bessel function [1] of order u {not to be confused with the 
symbol K we use for the kernel of V). 

In summary, the scalar kernels that we consider in this paper will always have the properties: 

(Kl) K is positive definite; 

(K2) K is symmetric, i.e. K{x) = K(—x), x G M.^ . 

In addition, in certain sections wc will introduce the following simplifying assumptions: 
(K3) K is twice continuously differentiable, K S C^(IR^); 

(K4) K is rotationally invariant, i.e. K{x) = 7(||a;||i{D), x G M"'', for some 7 G C^([0,oo)). 

Note that if (K4) holds then 7(0) > \-i{p)\ for all p > Q by (Kl) and (K2). Also, the bell-shaped 
Bessel kernels of the type (3) satisfy all of the above when A; > y + 1. 

Now fix any admissible Hilbert space of vector fields. The space i^([0, l\,V) is the set of func- 
tions f : [0, 1] — > V such that: 

/•I 1 

ibiii,p([o,i],v) •= (y )\\v < 00. 

The space i^([0, 1], V) is a subset of L^{[0, 1], V) and is in fact a Hilbert space with inner product 
{u, v) L2([o,i],y) := /q {u, v)v dt. It is well known from the theory of ordinary differential equations [7] 
that for any v G i^([0, 1], V), the Z)-dimensional non-autonomous dynamical system z = Vi{z). with 
initial condition z(io) = has a unique solution of the type z{t) ~ 'ip{t,tQ,x). Let <^^"j(.t) := 
^{t, s, x); fixing t = 1 and s = we get (p'" := ip^i, which is the diffeomorphism generated by v. For 
an admissible Hilbert space we will call the set 

Gv := {^^ :veL'{[0,l],V)} 

the group of diffeomorphisms generated by V; by [31, Chapter 12] it is a metric space and a topo- 
logical group. But, in the language of manifolds, Gv is not an infinite-dimensional Lie group [19]. 
V is not a Lie algebra, but is the completion of the Lie algebra of C°°-vector fields with compact 
support with respect to || \\v- 



2.2 Definition of the distance function 

For velocity vector fields v G L'^i[0, 1], V) and landmark trajectories q G Q define the energy 



E^[v,q] = E[v,q] := ^ (\\v{t, 



N 



|2 

\v ' 



^ll\\^it)-<t,q''{t)) 



a=l 



dt. 



(4) 



where A G (0, 00] is a fixed smoothing parameter (soon to be described). We claim that a dis- 
tance function d on {R^) between two landmark sets (or shapes) / = (a;^,x^, . . . ,a;^) and 
I' = {y^ ,y'^, . . . , y^) can be defined as 



d{I,r) := inf {/E[^ ■.veLmO,l],V),qeQ with q{0) = I, q{l) = /'}; 



(5) 
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in the next subsection we will argue that the above function is in fact a geodesic distance with respect 
to a Riemannian metric. We treat the minimization of (4) as our starting point; it is the "energy of 
a metamorphosis" as formulated in [31, Chapter 13]. 

The above infimum is computed over all differentiable landmark paths q € Q that satisfy the 
boundary conditions ((/"(O) = x"" and q°'{l) = y", a = 1, . . . , iV), and vector fields v G i^([0, 1], V). 
The resulting landmark trajectories {q'^{t),t G [0, l]}a=i,...,7v follow the minimizing velocity field 
more or less exactly, depending on the value of the smoothing parameter A G (0, cxd]; it is a weight 
between the first term, that measures the smoothness of the vector field that generates the diffeo- 
morphism, and the second term, that measures how closely the landmark trajectories actually follow 
the vector field. 

The exact matching problem is the following: given two sets of landmarks / = (a;^,a;^, . . . ,x^) 
and I' = {y^ ,y'^ , . . . , y^) with x"" ^ x^ and j/" ^ y^ for any a^h, minimize the energy 

^ooH := / \Ht, )\\ldt 
Jo 

among all v € L^([0, 1],^) such that ip'"{x°') = y"', a = 1, . . . ,N. In this case the landmark 
trajectories are defined as the solutions to the ordinary differential equations = v{t,q°'), a = 
1, . . . ,N. Note that this is equivalent to solving (4) for A = oo, since such equations are obtained by 
setting the integrands of the second term in the right-hand side of (4) equal to zero. When A < oo 
in (4) we have regularized matching, i.e. the landmark trajectories "almost" satisfy such set of 
ordinary differential equations; this allows for the time varying vector field to be smoother. For 
this reason the second term in the right-hand side of (4) is often referred to as sm,oothing term; 
by allowing smoother vector fields the distance d is made tolerant to small diffeomorphisms and 
therefore more robust to object variations due to noise in the data. 

2.3 Minimizing velocity fields and Riemannian formulation 

By manipulating expression (4) we will now show that it is equivalent to the energy of a path q € Q 
with respect to a Riemannian metric. 

Notation. Consider a landmark q = {q^ , . . . ,q^) in £^{R^). The D scalar components in Eu- 
clidean coordinates of the N landmark trajectories = ((7"^, . . . , g"^), a = 1, . . . , N can be ordered 

cither into an N x D matrix or in a tall concatenated column vector. We shall always use indices 
a,b,c, . . . G {1, . . . , A^} as landmark indices, and i,j, k, . . . € {1, . . . , D} as space coordinates in M^. 
We will associate to each of the A'' landmarks q'^ G R^^^ a momentum Pa gM.^^^ (defined in the 
next proposition) which we will write, in coordinates, as pa = (pai, ■ ■ ■ iPao), for each a = 1, . . . , TV. 
The components of momenta can also be ordered into an N x D matrix or in a long row vector. We 
chose superscript indices for landmark coordinates and subscript indices for momenta. 

For a given set of landmarks {q^,. . . ,q^) G we will define the symmetric N x N matrix 

K(g) :— {K{q°- — q'')) ^ j^- The matrix K(g) is positive definite by property (Kl) of the kernel. 

Proposition 1. For a fixed landmark path q = {g" : [0, 1] — )■ M'^j^j^ G Q there exists a unique 
minimizer with respect to v G ^^([0, 1], 1^) of the energy E[v,q], namely: 

N 

v*{t,x):=Y,Pait)K{x-r(t)), [0,1], a;eM^, (6) 

a=l 
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where the components of the momenta are given by: 

N 



PaAt) = E + • l^'^W, t G [0, 1], (7) 



6=1 



a = 1, . . . ,N, i = 1, . . . ,D (here Ijv indicates the N x N identity matrix). 

Remark. What the above proposition essentially says is that the vector field of minirniim energy 
that transports the A'' landmarks along fixed trajectories is, at any point of time, the linear combi- 
nation of N lumps of velocity, each centered at a landmark point. The directions and amplitudes of 
the summands are determined precisely by the momenta. 

Proof of Proposition 1. Using property (ii) of the admissible Hilbert space V , [31, Lemma 9.5] shows 
that for given q= (g^, . . . , g^) e £^(]R^) we have the orthogonal decomposition 

V ={v&V ■.v{q'')=Q,a = l,...N}(B{v = YJl=^aaK{ -g") : G M^}. (8) 
Thus the minimizer must have the form 

JV 

v{t,x) = Y,aa{t)K{x-(f{t)), t e [0,1], a;GM^, (9) 

a=l 

for some coefficients aa G C([0, 1],M^), a = 1, . . . , A'', to be computed. For velocities of the type (9) 
the energy (4) can be rewritten as 

I D N 

E[v,q] = - i>)au + A \(XaiK{it - f) - f dt. (10) 

•^0 i=l a,b=l 

Setting the first variation of (10) with respect to coefficients aai to zero yields the momenta (7). □ 
It is convenient, at this point, to introduce the ND x ND, block-diagonal matrix 
/(K(9) + ^)"' ••• \ 



9{q) ■■= 



V ••• (K(g) + ^) 7 



(11) 



where the N x N block (K(g) + x) ^ ^® repeated D times; the choice of symbol g is justified by 
the fact that (11) is, as we shall see soon, precisely the Riemannian metric tensor with which we are 
endowing the manifold of landmarks, written in coordinates. 

Thus for a fixed path g G Q the minimizer of E[v,q] with respect to v £ -^^([0, 1],!^) is given 
by (6); since it depends on q we will write it, with an abuse of notation, as v*{q). We can define 

E[q]~E[v*{q),q], (12) 

which depends only on the arbitrary path q & Q. The energy (12) is "equivalent" to the en- 
ergy E[v,q], in that: 

(a) if (w,g) minimizes E[v,q] then q minimizes E[q], and E[v,q] = E[q]; 

(b) if q minimizes E[q\ then {v*{q),q) minimizes E[v,q], and E[v*{q),q] = E[q]. 
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Proposition 2. For an arbitrary landmark trajectory q G Q the energy E[q] is given by: 

m= f'q{tf9{q{t))mdt= f (t) ci>\t)U{q{t)) + \y' dt (13) 

In the above equation q{t) is intended as an TVD-dimensional column vector obtained by stack- 
ing the column vectors (g^*(i), . . . , g^*(i))-^, i = 1,...,D (again, the superscript ^ indicates the 
transpose of a vector). 

Proof. Following definition (12), formulae (7) for the momenta are inserted into the modified expres- 
sion (10) for energy E[v, q]. Simple matrix manipulations finally yield the right-hand side of (13). □ 

Remarks. Expression (13) has exactly the form of the energy of a path q with respect to Ricmannian 
metric tensor (11). Whence given two landmark configurations / and /' in >C^(M^) we have that 
if q minimizes (13) among all paths in q € Q such that q{0) = I and q{l) — I' then {E[q]y^^ is the 
geodesic distance between / and /'. By point (b) above we also have that {v*{q),q) is a minimum 
of energy E[v,q], so d{I,I') defined in (5) coincides with (£[g])^/^ and is the geodesic distance 
between I and /' with respect to the metric tensor g. 

The Lagrangian function that corresponds to the energy (13) is: 

C{q,q) = -q^g{q)q=-^ ^ E '^"^'^ (^l?) + ' ^^^^ 

o,6=l i=l 

In Hamiltonian mechanics [3, p. 60] the "momenta" are defined as Pai = dC/dq'^^, or, in vector 
notation, = dC/dcf^'^ (for i = 1,...,D). Applying such definition to (14) yields precisely 
equations (7) of Proposition 1. Whence the use of the term momenta is justified. 

Note that for small values of the parameter A the metric tensor g, written in coordinates, gets 
close (up to a multiplicative constant) to the ND x ND identity matrix; in other words, for A ^ 0, 
g converges to a Euclidean metric and the geodesic curves become straight lines. On the other 
hand, for A ^ oo (exact matching) the metric converges to [diag{K((3'), . . . ,K((7)}]~^ (block 'K{q) 
is repeated D times). In general, the block-diagonal form of the metric tensor g given by (11) 
follows from the fact that the operator L in (2) is applied separately to each of the components of 
the velocity field; however the dynamics of the D dimensions of q are not decoupled since all ND 
components of q appear in each diagonal block of g. 

In the case of exact matching landmarks "never collide" (their trajectories are precisely defined by 
diffeomorphisms of M^): it takes an infinite amount of energy to make any two landmarks coincide. 
So under the condition A = oo the manifold of landmarks can actually be taken as the set: 

C^{M.^) = {(P\ . . . , P^) I G M^, P" ^P^ lia^hy (15) 

Figure 1 shows the qualitative behavior of geodesies in £^(]R^), with A = oo. In the case 
illustrated on the left-hand side both landmarks travel in the same direction (from left to right, as 
indicated by the arrows): the two arcs of the geodesic "attract" each other, or in other words the 
two landmarks tend to "carpool" by using a velocity field with the smallest possible support so to 
minimize the L'^ part (i.e. the first term) of the Sobolev norm (2) of the velocity field. On the other 
hand when the two landmarks travel in opposite directions (as illustrated on the right-hand side of 
Figure 1) they try to avoid each other so that the higher order terms of the Sobolev norm are kept 
small; we shall return on the issue of obstacle avoidance at the end of this paper. A typical geodesic 
in £^(M^) (again with A = oo) is shown in Figure 2. 
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Figure 1: Two trajectories in Bullets (•) and circles (o) are the initial and final sets of 

landmarks, respectively. The grids represents the two corresponding difFeomorphisms ipQi. 

Conclusion. We have shown that distance d{I,I'), 1,1' E C^{M.^) defined in (5) is in fact the 
geodesic distance with respect to a Riemannian metric. In coordinates, the corresponding Rieman- 
nian metric tensor is given by (11), which is such that each element of its inverse (the cometric) 
depends on at most 2D of the ND coordinates. Whence the first and second partial derivatives of 
the cometric have a very sparse structure. This gives us motivation for deriving a general formula 
for computing sectional curvature in terms of the cometric and its derivatives in lieu of the metric 
and its derivatives, which will be done in the next section. 



3 Sectional Curvature in terms of the Cometric 
3.1 Generalities and notation on sectional curvature 

Let Ai be an n-dimensional Riemannian manifold. If we consider a local chart (U, ip) on the manifold 
with coordinates {x^, . . . , a;"), we have the induced 1-forms dx^, . . . , dx" and coordinate vector fields 
{di :— gp:, . . . ,dn — gfrr}- The metric tensor g : TM. Xyvi TM — > M can be represented as g\u = 
g{di, dj) dx' (J5 dx^ =: gij dx^ dx^ (here, as in the rest of the current section, we are using Einstein's 
summation convention). For each p £ M we get a positive definite matrix with elements gijip) = 
gp{di, dj). With an abuse of notation we will write gij{x) instead of {gij o tp^^)[x), x G f{U) C M". 

Notation. We shall denote the partial derivatives of the elements of the metric tensor g as gij.k{x) := 
^gij{x) = dkgij and gij,ki{x) ■= gJl^k 9ij{x) = did^g.^j. Also, we wiU indicate the cometric 
g~'^\u — g^^di (8) dj (so that g'-'gjk = S^) and their partial derivatives with g^-' j^{x) := gffc5*''(2;) 
and g"\kei^) g'^x)- 
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Figure 2: A typical geodesic trajectory in Bullets (•) and circles (o) are the initial and final 

sets of landmarks, respectively. The grid represents the corresponding diffeomorphism (p^i. 

For a tangent vectors X = X^di we consider the l-form X^ := X^gijdx^ —: Xjdx^ (indices 
lowered), and for a l-form a = aidx^ we have the tangent vector a' := aig^^dj (indices lifted). 

Indicating with X{M) the space of smooth vector fields on the manifold A^, let V : X{A4) x 
X{A4) X{M.) be the Levi-Civita connection [16, 21] of the Riemannian manifold. The Christof- 
fel symbols T^^ are defined by Va^c^j = T^^dk^ and it is well known that they have the form: 
r^j- = \g^^{gu.j + gji,i ~ gij.e)- The Riemannian curvature endomorphism is the map R : X{A4) x 
X{M) X X{M) X{M) given by R{X,Y)Z = ^x^yZ - VyV^^ - "^\x.y\Z. In local coordi- 
nates R{di,dj)dk = Rlji^di, and Rijkm '■— {R{di,dj)dk,dm) g = gmiRljk- The Riemannian curvature 
tensor acts on vector fields as follows: 



R{X, y, Z, W) := {R{X, Y)Z, W)g (16) 

and in coordinates it is written as i? = Rijkmdx^ (E) dx^ (E) dx^ (E) dx™. The Riemannian curvature 
tensor has a number of symmetries: (i) RijM = -Rjike] (ii) Rijki = -Rijek] (hi) Rijki = RkUj] 
and (iv) Rijki + Rjku + Rkijt = (first Bianchi identity). With such conventions, the sectional 
curvature associated to a pair of non-parallel tangent vectors X and Y is computed by: 

ir(XY\ = RiX,Y,Y,X) ^ R^.kmX^Y^Y'^X"^ 
^ ' ' \\X\\l\\Y\\l-{X,Y)l \\X\\l\\Y\\l-{X,Y)l ^ > 



In order to express the numerator of sectional curvature (17) in terms of the elements of the 

fc , and g'\^^) 



cometric and its derivatives (i.e. g''^ , g^'' k^ and g*"* j,^) we consider the covariant expression of the 



Riemannian curvature tensor: 

H .— Ktjkm g g g g , (.loj 
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which we call the dual Riemannian curvature tensor. Similarly we consider the covariant or dual 
Christoffel symbols 

■■= rg'^QkuTlj , (19) 

which are symmetric in the indices r and s. 

To achieve notational compactness we will use the following symbols: 

g'''' ■■= g'\i 9^' and g'''''':=g'\^^g^'g^'; (20) 

Using that g = implies dkg = —Q~^ ■ d^Q ■ one immediately sees that 

= - IguM""'" + 9^'"'' - Q'^'"')- 

Proposition 3. The following expression holds for the Riemannian curvature tensor: 

'^Rijkm — 9ik,jm H" ^jm^ik 9jk,im 9irnjk H" ik^ jmdoe^ jk^ im9oi./3 ' (^1) 

For a proof see [24, §24.9]. 
3.2 Mcirio's formula 

Proposition 4. The following expression holds for the dual Riemannian curvature tensor: 

(2j^ursv gUS,rv grv,us _|_ grs,uv _|_ gUV,rs _|_ t2Y'rvY'US gpa 2P'^^p^^^/^^ 

I „rA,ti „IJ.i>,s „r\,u „t^s,v , _tiA,r nlJ.s,v „uX,r „liv,s fOO\ 

+ 9 9\p. g g .9a^ g'^ + g gxp. g - g g\ng (^^) 
+ g g^i^g + g g\t^g ~ g g\ng - g gxfj.g ■ 

Proof. We will manipulate (21) and write it in the form Rijkm — giugjrgksgmv 

^urev factor- 
ing giugjrgksgmv out of cach term; what will be left will be precisely the expression for i?"'"*"'. 
The terms in (21) involving Christoffel symbols are, by (19): 

Tf,r^M = 5iu5fe.5"" rr gjrgmvg'^ g^^ = giugjrgksgmv (r^rr^^'^) , (23) 

and similarly: ^jk^Lg^ff = giugjrgksgmv (r-r^g^-) . (24) 

As we noted before, if 5 = then djg = —Q~^ ■ djQ ■ and similarly it is the case that 
9m9jg = ■ {dmQ ■ ■ djQ + djQ ■ Q~^ ■ d^Q - d^djO) ■ Q~^, i.e., in index notation, 

gikjm — giu[g ,m.g\n9 ,j +g jgxng ,m -g ,jm)gsk 
= giugksS|^l{g^\gx^.g''\^ +5"^^5A^5''% -5"^^^) 
= giugksgjrgmv[g'^^g''^{g''\r, gx^g^^i + 9^\i 9x^g^',r> -9^', in)] 



= giugjrgksgmv [g g\n g^ + g gx^ 5 - g j > (20) 

where we have used definitions (20). Similarly, we can achieve the factorizations: 

gjm,ik = giugjrgksgmv {g'^^^^gx^ 5""'^ + 5''^'^5a^ s""'" - 5™'"^ , (26) 

-gjk,im = giugjrgksgmv ( - 5"^'"5a^ .g"^'" - 5"^'"5am 5'^'^'" + /''"I , (27) 

-gim^jk = giugjrgksgmv ( - ff"^'^5A^ 5""'^ - 5"^'^5a^ g^""'' + 5""'") • (28) 



Inserting (23)^(28) into (21) we can write Rijkm = giugjrgksgmvR'^'^^'" , with i?"™ given by (22). □ 
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Proposition 5. The dual Riemannian curvature tensor may also be written as follows: 



in — — g ' — 9 -r 9 



^us^rv ^rv^us \ „rs,uv _|_ gUv^rs 



+ ^{5™,p -5™,p (5''"'^ +9"^'") (5"''''^ +5""''')} 



(Ti) 

(T2) 

(T3) 
(T4) 
(T5) 

(Te) 



Proof. We will expand and recombine the terms in expression (22). The terms involving second 
derivatives need no manipulation and correspond to term Ti. The terms in the second line of (22) 

can be written as: 

9 9\ii9 -9 9\fj.9 +9 9\ij.9 9 9\n9 —{9 9 )9\n{9 9 ) 



which is precisely Tg. It is also the ease that 



^^p 9 -9 9\f^9 " 



2ia ,p; 

Similarly one can prove that: ~2r„^T'^" gP"^ + g^^^'gx^ g''"''" + 5"^'"5Ap 5^"^'' = Ta + T4. 



T3 + Ks^'-'" - 5^''''^)5A^(5''"'^ - 3'^^'") = T3 + T5. 



□ 



For any point p € M. and an arbitrary pair of tangent vectors X = X^di, Y = Y^di in TpM. we 
consider the covectors X'' = Xidx^ and F'' = Fjda;' in Tp with Xi = 9ijX^ and Yi = gijY^ . The 
numerator of sectional curvature (17) may be rewritten as Rij^rnX'^Y^Y'^X"^ = R^^^^ X^Y^YsXy. 



Theorem (Mario's formula). For an arbitrary pair of vectors X = X'^di and Y = Y^di in TpA4 

the numerator of sectional curvature (17) at point p G Ai may be written as: 



g{RiX,Y)Y,X) = 7?"'-^"X„y,y,X„ = 






= ix„Yr - Y„x,.) (i'/'"-'"'' + 




i'/"-''flA„fl""'')(Xl^-KX„). 



Moreover, if we extend X'' and Y^ locally on M. to constant 1-forms in terms of local coordinates 
(i.e. make its coefficients Xu,Yr constant functions), then the formula becomes: 



g{R{X,Y)Y,X) = 

= [^xxiWY^f) + iyy(||x^||2) - i(XF + FX)ff-l(x^y^)} 

+ {l\\d{g-\x\Y'))r - y-'{d{\[x'r),d{\[Y'r))} - U[X,Y], [X,Y]), 



11 



where the term, in the first set of braces equals the sum of the first two term,s in the coordinate form, 
the term in the second set of braces equals the third term in the coordinate form and finally the last 
terms are equal. In the above formula, H-'^^lP = X^X^^™ and = YrY^g^^. 

Proof. We will write the six terms provided by Proposition 5 as T"''*^, i = 1, . . . ,6. We have: 

rpursv V V V V us,rv V V V V rv,us y v v y i ,,rs,uv y v v y i ,,uv,rs y v v y 

_ g^^'^"" (^ — X^ijY^YgXy — X^YyYyXg + X ^YyY g X y + XyYjYyXg) 
_ ^"«>™(X„1^ — YuXr.){XsYy — YgXy), 

where the second step follows from relabeling the indices. As far as T2 and T3 are concerned, 

(Tr-+T«'^-)X„y,F,X, = -\{Y,Ygg-\p gP- g^\„X^X, - Y,X,g^\^ gP-g^^X^F,} 

+ \{YrYsg''\^[gP^^^ ^gP^^^)X^X, + X^X,g^\^[gP-^' +g'''-'-)YrYg 
- YrX, g'^'.pig'"'' + g'"'")X^Ys - X^Y, 5"^,^ (s"'^''' + gP^''')Yr.X,} 
= - l{2YrYgg'-\p gP-g^\,X^X, - 2YrX,f\p gP^ g^\^X^Y,] 

+ \ { 2YrYg g^\^ g'^'^'^X^X, + 2X^X, g^\, g'^'^'^YrYg - 2X^Yg g^^^ {g^^^- + g"^^^) Y^X, } 

= — \ 9^^ g'"^ 9^^ ^a{YuYsXj.Xy + YrYyX^Xs — Y^XyX^Yg — YyXrXgY^^ 

+ ,P 9'^^'^ {^uYgXrXy + XuXgYj-Yy — XyYgYj-Xy — XgYyVyX^} 

= ( - 3 .9"';. ,9™'" + .9"^p 9"'''") {XuYr - YMiXsYy - YgX,), 

where, once again, step (*) follows from relabeling the indices. Also, one can easily see that 
Tur^vx^YrYgY, = -lYrYgig^^'^ - g^''')gxM'''" ' = 0. Finally, 

= \ YMg^^'- - g^^'ngxM^-' - 5^''")^un + YrX^{g^-^^ - g'-^ngxM"-' " 

= I YrXu{9^^'^ - g^'^ngxM'''' - g'^'nXvYg 

= I YrX^g^^^^gx^g'^-'^ - g^^'^^gx^g^'^^ - ^^"'''flA^s''"'^ + 5^"'''5A^5''^'''}^.n 

= — I g^^'^ g\lJ,g^^'^ { — Y^X^XgYy + Y^X^XyYg + Y^Xy^XgYy — YrXuXyYg^ 

= - i 5^"'''5A^5''^'''(^„n - Y^Xr){XgYy - YgXy). 

Divide by 2 to get the coordinate formula. The non-local version of the formula follows easily by 
bringing the X and Y's into the formula. Thus (indicating di with the subscript 

y„x,(/"'™ + = XrXy{YgY^ g'^g^^ + ysYu g''',p g'^a 9"") 

= XrXy{{\\Y^f),p,gP'g-'' + {\\Y^ f) ,p g'^, 5"") (because Yg,Y^ are constants) 

= Xyg'^^{Xrg''^{\\YY),p)^^ = X-(X^(||F^f ),,)^^ = XX{\\Y'f). 

A typical term from the third part of Mario's formula is rewritten like this: 
Y^Xrg'^\,g'^^'^YgXy = XrXy{\\Y'r)^^g^^'- = {\\Y'r)J\^^^ 
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the other terms are similar. Finally, it is the case that: 

= {{Xug^^),r,Y^ - {Y^g^^),r,X^)dx = {x\Y^ - Y\X^)dx = -[X,Y], 



and the proof is easily completed. □ 
Remark. It is convenient to split Mario's formula in four terms: 

Ri := liX^Yr - y„X,)5««.™(X,K„ - Y^X,), (30) 

i?2 := U^uYr - Y^Xr)g''\p gP^-'"{X,Y, - (31) 

i?3 := \{X^Yr - YuXr) ( - \ .9"% .9™'") {X^Y, - Y,X,), (32) 

:= i {X^Y, - Y^X,) ( - I g^^^-gx^. <?''^'") {X,Y, - Y,X,) ; (33) 



all the terms with the exception of i?4 (where g appears, but not its derivatives) depend only on 

elements of the comctric and their derivatives. 

Remark. The denominator of sectional curvature (17) can also be expressed in terms of the comet- 
ric: 

ll^llgll^ll^ - = - g'^'^gn- (34) 

4 Curvature of the Manifolds of Landmarks 

In this section we will apply Mario's formula to the computation of sectional curvature for the 
Riemannian manifold of landmarks, introduced in section 2. We first introduce the Hamiltonian 
formalism, since it will allow us to write the geodesic equations in a simple form and to introduce 
geometric quantities that will eventually appear in the formula for sectional curvature. 

4.1 Hamiltonian formalism 

On the A''Z)-dimensional manifold £ = {MP) of landmarks wc consider the Riemannian metric g 
given, in coordinates, by the matrix (11); it is in block-diagonal form and we write its generic 
element as g(ai){bj)i with n, 6 = 1, ... ,7V (landmark labels) and i,j = 1, . . . ,D (coordinate labels, 
respectively of landmarks a and b). More precisely: the matrix g{q) is made of D square {N x N) 
blocks; indices i,j = 1, . . . ,D indicate the block, whereas indices a,b= 1,. . . ,N locate the element 
within the (i, j)-block. Therefore if we indicate with habio) the generic element of the N xN matrix 

(K(g) + ^)"^ we have that 

g(ai){bj) = Kb{q)Sij, a,b=l,...,N, i,j = 1,...,D, 

where Sij is Kronecker's delta. Similarly, if we indicate as gi^"') the elements of the cometric 
tensor S'(9)"S they are given by 5r(")(^-?)(g) = h''^{q)S'^, where h"^{q) = K{q°'-q^) + ^. In analogy 
with the notation introduced in section 3 we also denote the partial derivatives by g^"'^^^''^^ ^^^^ = 
^g(-i)m and g^'''^^''\^,k)(de) = they will be computed later. 
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For simplicity from now on we shall assume that A = oo . i.e. that wc are dealing with exact 
matching of landmarks so that (M.^) has the form (15). The element of the cometric becomes 
gr(a»)(6j)(g,) = K{q°' — q'') 6''^ and the Hamiltonian [16, p. 50] for the system can be written as: 



N D , N D 



nip,q) = lp^g{q)-'p=lj2 J29^'''^^''\'l)PaiPl>i = lj2 J2K{q'^-q')6'^p,,pi,j 



a,6=l i,i = l a,6=lj,i=l 



that is 'H{p,q) = \Y. ^(9" ' <l^){Pa,Pb)^D- 

a,6=l 

Proposition 6. Hamilton's equations for the Riemannian manifold of landmarks are: 

N 

a=l,...,N. (35) 



Pa = -J2'^K{q''-q''){pa,Pb)^ 



6=1 



Proof. Equation (7) can be written as q"^ = J2b=i I^i^f ~ Q^)Pbi, for a = 1,. .. ,N, i = 1,.. .,D; 
alternatively, computing g"' = ^7 yields the same result. Also: 

^K{q^' - q'\ . . . , g*"^ - g'^^) = Ef=i i^(9' - ' l"') 

so that 



Ef=i - Q") {^l -m = M-^ - 1") (Si - K) (36) 



Pai = -^{p,Cl) = -i Ef=l If (9" - 9^) {Pa.Pc)^o + i Y.tl Bil' - I'') iPb,Paho 

= -ELil^{<l''-Q''){Pa,Pb)uo; 
in (*) we used the skew-symmetry of VK{q°' — q^) in indices a and h, which follows from (K2). □ 
Corollary 7. IfPa{to) = for some landmark a = 1, . . . ,N and time to € K, then Pa{t) = 0. 

4.2 Notation 

Prom now on we shall also assume that (K3) holds, i.e. that the kernel K is twice continuously 
differentiable; for the time being we will not assume rotational invariance. We define: 

K"'' := Kiq" - q'') G R, 

8K 

diK{x) —^(x), diK'^'' := diK{q'' - q'') G M, 



:= {diK, ■■■ , doKy , VK"" := VKiq" - q") e M^, (37) 



D^'K := Hessian(ii'), D^'K"" := D^K{<f - q") e 
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Note that VK'''' = -VK'"', VK"" = and D^K"'' = D^K^"", for all a, 6 = 1, . . . , iV, by (K2). 

For a fixed set of landmark points g in £ = £^(K^) consider any pair of cotangent vectors 
a,l3 G T*C: we shall write a = (ai, . . . ,q:jv) and /3 = . . . ,I3n), where each component is D- 
dimensional. We define the vector field 0;'^°'' : MP — )• and its values at the landmark points by: 

N 

6=1 

which are, by virtue of formula (6), the velocity field a^"'' on induced by the landmark mo- 
mentum a = (ai,...,ajv) and the corresponding landmark velocity a" e TqC (which obviously 
coincides with the first of Hamilton's equations (35)). Note that a' = (aj, . . . ,q5v) is the tangent 
vector in TqC with metrically lifted indices. Note that a^°'' is the horizontal lift [10, p. 148] of the 
tangent vector a" on the admissible Hilbert space V: simply put, of all vector fields v : — )• 
in V such that v{q°') = (a")", a = 1, . . . , N , a^"" is the one of minimum norm. 

The curvature of the Riemannian manifold of landmarks will be expressed in terms of three 
auxiliary quantities which we now introduce. We will call these force, discrete strain and landmark 
derivative. We start with the force. For a fixed covector a = (ai, . . . ,ajv) € T*£, having the dual 
vector extended to a vector field a^""" on all of allows us to take its derivatives at the landmark 
points, a D X D matrix- valued function on M^: 

N 

6=1 

N 

(Z)a^-)j(g») = ^a,/^«V-. 
6=1 

For a trajectory {q{t),p{t)) of the cotangent flow one has that {pi{t), . . . ,pN{t)) G '^q(t)^ ^'^^ * 
where the trajectory is defined, so the above notation can be used to rewrite Hamilton's equations 
in a more compact form. In particular, the following result holds. 

Proposition 8. The second of Hamilton's equations (35) can be written as 

Pa = -Dp''°''{q'')-Pa a = l,...,N. (38) 

Proof Pai = - J:tl^^K''''iPb,Pa)Rn = " Ef=l {E^l ^^K^'' Pb,)Paj = " Ef=l Pa, = 

-(Z)p^°"(g") -Pa)., for anya = l,...,7Vandz = l,...,£). □ 

For a fixed cotangent vector a € T*£, this motivates defining the negative right-hand side of (38) 
to be force: 

Fa{a, a) := i?a^°'^(g") • a„, a = l,...,N. 



N 

a^°^{x) := 

b=l 

{J)" := a^°'{q"-) = 
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The full bilinear, symmetrized force may be thought of as a map F : T*C x T*C — >■ T*£. We call 
the covectors given by this the mixed force, with the definition: 

F,{a,l3) := l{Da^°'' {q'^) ■ /3„ + • a„), 

j=l 6=1 6=1 

for a = 1, . . . , A?' and i = 1, . . . ,D. (The angle brackets are inner products in M^.) Note that the 

"complete" cotangent vectors a = {ai, . . . , a^v) and /3 = (/3i, . . . , Pn) (not only their a-components) 
are needed to compute each component Fa{a,f3) of the mixed force. The mixed force has simple 
interpretation. If we extend a and /3 to constant 1-forms on £, then the differential of the map 
q ^ g-^{a,(3) = J2a,b^il°' - Q''){(^a, l}b)Rr, is given by: 

N D 

a, 6=1 2 = 1 
N D 

= J2 - ((«a'/3i.)K° + (^a,"^)^^) = 2F{a,l3). (40) 

a, 6=1 2=1 

For a fixed a &T*£ we define the discrete vector strain: 

N D N 

S^^ia) := (a")" - (a'')^ or S"*^(a)* := ^ Y^iK"" - K^'')5'^ a^j = Y^iK""" - K^'')aci 

C=l J = l c=l 

for all a,b = 1,. . . ,N (we call it like that because it measures the infinitesimal change of relative 
position of the landmarks a and b induced by the cotangent vector a). These are vectors and are 
skew-symmetric in the points a,b: S"'^{a) = —S^"'{a), S°'°'{a) = 0. The scalar quantities: 

N D 

c=l i=l 

we define to be the scalar compressions felt by kernel K; they are symmetric (since both factors in 
the inner product are skew-symmetric), i.e. ^^{a) = C^°'{a), with the property C""(a) = 0. We 
call these compressions because if X is a monotone decreasing function of the distance from the 
origin (the most common case), then V-ftT"'' points from to q''. 

Finally, if v and w are any two vector fields on the manifold of landmarks, we may write their 
Lie derivative as the difference of covariant derivatives: 

b,«;]z: = Vp-*(ti;)-V^'«-*(^) 

where the flat connection on C is just the one induced by its embedding in M^^. In other words, 
Y^flat^^-^ is the usual derivative of w in the direction v if we use the coordinates on landmark 
space: that is, V^^''^''*(w;) := Z^^j v(w'"*)9aj = J^ai^bj '"''^ i^bjw"'')dai- If a, P are constant 1-forms 
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everywhere on we can take w = a' and w = ji'^, now as vector fields on C, and then we find: 

a,i b,j a,z b,j c 

= E E ^j-^" - St) Pa 9« = E E - 9« 

a.z b,c.j a.i b,j 

= E E - v^"''>M«/3fi = E (E 

This is a vector in TqC which we define to be the landmark derivative of /?" with respect to a". The 
coefficients with respect to dai , daD (for fixed a) are the elements of the following vector: 

N N 

D-{a, p) := E C"'(«)/36 = E ' ^'')("c, '^K-^)^op^ , a = l,...,N. (41) 

6=1 6,c=l 

We have that D{a, p) = (r>«(a, is the TVZJ-dimensional vector of the coefficients of Vf J*'''*(^'*) 

with respect to the basis of TqC. In particular, the coefficients of the Lie bracket of a' and /3' 
as vector fields on C are given by -D(a, P) — D{fi, a). 

4.3 General formula for the sectional curvature of C^{EP) 

Wc can write sectional curvature of C'^ {^P) in the following way, where we have split it in the terms 
introduced by (30)-(33). 

Notation: from now on ( , ) will indicate the dot product in 'MP , while { , )tc and ( , )t-c 
will be the inner products in the tangent and cotangent bundles of £ = £^(M^), respectively. 

Theorem 9. The numerator of sectional curvature of {R^), for an arbitrary pair of cotangent 
vectors a and p, is given by R{a^ , p^ , P^ , a^) = X]i=i with: 

^1 = 5 E ® ^"''('^) - ^« ® 5<^*(a))^(lD D'^K''^) («(, ® 5«*(/3) - pt,(S) 5"^(a)), (42) 
^2 = E ((^"(a,a),i^a(/3,/3)) + {D\fi,l3),F,{a,oc)) - {D\a,li) + D\p,a),F,{a,f5))), (43) 

a 

i?3 = \\F{a,f3)\\l,^ - {Fia,a),F{p,0))^,^ 

= E P)^Fc{a, (i)) - {Faia, a),F,{P, /?))) , (44) 

ac 

i?4 = -l\\[aKp%\\l^ = -l\\D{a,p) - D{p,a)f^_,. (45) 

In the formula we have used the definition: (t;i(g)t;2)^(Mi(g)M2)(u!i(S)t«2) := {vf Mi'Wi){v2 M2W2) 
for the first term while we have used the norm for DxN matrices || J||^ := X^^i ^^^6=1 JiaJib-^ab 
for the fourth term R4. 

The theorem is proven by applying Mario's formula to the cometric of the manifolds of landmarks. 
One needs to compute the elements of the cometric and its derivatives in terms of the kernel and 
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its derivatives (37). In agreement with notation (20) we will define (note that we will keep using 
Einstein's summation convention wherever possible): 

{ai){bj),{de) (ai)(bj) {ck){M) J (ai){b3),{ck){de) (ai)(6j) {ij.p){ck) {ia){M) 

^ ' ^ tick) y y • y ,ii^p){$<^) ^ ^ 

Lemma 10. It is the case that 

.9^"'*'"l(cfe) =5fcif"^<5»-<5,^)<5^-'', (46) 
9^'''^^''\ick)W = duK"' (S^. - S'^J (S^, - 5^) (47) 

g{ai)(bj),{M) ^ Q^j^ab ^j^ad _ j^bd-^ g^J ^ ^^g-j 
g{ai){bj),{ck)idi) _ Q'^^X"'^ {K"" — K'"^) {K°''^ — K'"'') 6^^ . (49) 

Proof. Since = K'''>S'^ and also ^K{q'' - q^) = dkK''^{S^ - d^) by (36), equation (46) 

follows immediately. Similarly to (36) one can prove that g^dkK{q'^ — q^) = dJ^K^^ {5'^ — 5^), 

whence: = fl^"^^''\(e.) = dfkK''' i^d ' S'dHS^c ' ^l) so (47) holds too. Now, 

by expression (46): 

g(ai){bj),{dl) ^ gick){de) ^ Q^j^ab (^a _ ^6) ^ij j^cd §ke ^ Q^j^ab (^j^ad _ j^bd^ §ij _ 

which is (48). We can use (47) to compute g(ai)m,{ck)(de) ^ ff^"*^^*"^^ (^^j(j^) g(tip){ck) g(ia){di) . 
g{ai){bj),{ck)ide) ^ Y.f,p^a dpaK"'' {6^ - S''^) ((5| - <5|) S'^ K'"' dP'' Ki<^ 5"^ 

which completes the proof. □ 

Proof of Theorem 9. We will compute terms -Ri , . . . , i?4 introduced by formulae (30)-(33) . For sim- 
plicity, sometimes wc will write Dajj™ instead of Da^°^ {q""). 

• Computation of Ri. We have i?i = \{aaulicr- PauCtcr) g''°''^''^^^^''''''^^'''^'"\oi},slidv- PbsQ-dv)- Inserting 
expression (49) into such formula yields: 

2^1 = Eall indices ("au/3cr - Pauacr) d^^K'^'' (K'^^ - K'^'^ (^"'' " K''^) <5"' {absPdv " Pbsadv) ■ 

Performing the above multiplications gives rise to four terms, which we will now compute one by 
one. First of all we have: 

2^il,l := Eall indices C^auPcrabs^dv d^^K^^ (R-^ - K^^) (if - K^'') 5^^ 

= Habrv [T.us0^au6-'abs] [Ec(^"' " K"^)Pcr\ d'^^K^' [E<i(i^"' K''^)Pd.] 

= Eab oZo^b E™ S-\pr dj^K-^ S-\pr = Eab {S'^\P)fD^K-'' S-\/3); 
= E„6 ("a ® S'''{I3)) ^ (Id ® D^K-") {aj, ® 5''^(/3)) , 
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where, once again, the superscript indicates the transpose of a vector; similarly, 

2i?i,2 := - Eaii aauPcrhsadv d^^K'^'' (K'^' - K''^) (if "'^ - K'"') 6^' 

2i?l,4 := Eall PauQ.crPbsadv (if " K^^ (if - if 

Now we can take the summation i?i = Ei=i which yields precisely expression (42). 
• Computation of R2. We may combine equations (46) and (48) from Lemma 10 to get: 

f ^""■'''"Vp) i?^^''^^''''^'^'*''^ = Eap dpK"'' {SI - (5*) (5"" a„if ^= (if - if'='') SP^ 

= drK"^ [d^K'"'{K'"^ - K""^) - dyK^'={K^'^ - K'='')] S'^". (50) 
Inserting (50) into 2i?2 = (aauPcr - /3a«acr) 9^^''^^""^'^''"\c,bsPdv - Pbs^dv) yields: 

2^2 = Eall indices { ^auPcrabsfidv OrK^" [d^K^^K^" - K^") - d^K'^K"" - K^')] <5«^ 

- Pau^MPdv drK-'[d,K-'{K-^ - if^^) - d^K^'^K''^ - if^^)]<5"^ 
+ (3auaer(3bsadv drK^' [d.K^^K'^'' - K-'') - d,K'>-{K'"' - if^'')](5"^}, 
which immediately implies: 



i?2 = 



\ Eabcdi^a, at) {00, Vif«*) [{Pd, VK-%K»<^ 




{0d,VK^''){K'"^- 


-if 


{.=■■ 


R2,l) 


-h^abcdi^a, mPc, Vif«*) [{ad, Vif-)(if«'* 




{ad,WK^%K^'^- 


-if 


{=■■ 


i?2,2) 


- \ Eabcd (Pa , «b) (ae, Vif [{^d, Vif^) {K^^ 




(/3d,Vif^^)(if^'* 




{=■■ 


i?2,3) 


H Eabcdil^a, lib){ac, Vif"") [{ad, Vif-)(if 




(ad,Vif'"=)(if'"^ 




{.=■■ 


i?2,4) 



We will now manipulate terms ii2,ir • • 1^2,4 one by one. Since Vif" = — Vif , by relabeling the 

indices wc have 

^2,1 = Ea6cd(«a,a6)(/?c, Vif»^)(/3d, Vif-)(if«'* - if^'') 

= Ea6c("a, ab){Po, Vif (Ed^"'/3<Z " Ed^^'^<i, Vif^) 

= Ea6c("a'«6)(/3c, Vif«'')(5-(/3), Vif-) = E„6c("a' "6)(/3c, Vif «'')C-(/3) 
= Eabi^a, «6)(Ec^?"'(^)^c, Vif«*) = ^abi^a^ ",)(i?«(/3, /3), Vif«^) 

= Ea6 prVK-t>alaa = Ea D-{p, pYDal-^ ■ a„ = Ea(^"(/3. P), Fa{a, a)). 
Similarly, i?2,4 = '^^{D"{a,a), Fa{(i, (3)) ■ It is also the case that 

R2,2 i Ea5c("a, /3b){l3c, Vif »^) [(Ed (^"' " , ViC^^) - (E<i(if '"'-i^^'')^^, Vif^)] 

= -5 EabMa, Pb){l3c, Vif »^) [(5-(a), ViC"^) - {S^-{a), Vif 
= - 5 E„6c("a, Pbm, VK-"} [C-(a) - C''<=(a)] ; 
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relabeling the indices (and using the fact that VK""^ = —VK^"") yields: 
R2,2 = - \ Eabc Pb) + {ab, /3a)] Wc, VK-^)C-%a) 

= - 5 Eafc Pb) + {ab, Pa)] (EcC"''(")/3c, Vi^«^) 

= - 5 Eafa [{aa, Pb) + {ab, Pa)] {D^ia, p),VK-^) 
= - \ ZabD^ia, [VK-'p^aa + VK-'aJpa] 

= - \ EaD'^ia, Pf [Dp^°- ■ aa + Da'^ ■ Pa] = -Ea(^"(«> P),Fa{a, p)). 
Similarly, i?2,3 = -Y,^{D''{P,a),Fa{P,a)) . By the symmetry olFa{-,-), 

R2,2 + i?2,3 = -Ea(^"("> /3) + D-iP, a), Fa{a, P)) 

Adding the above sum to the expressions for ii2,i and i?2,4 finally yields (43). 
• Computation of R3 . We have 



But by Lemma 10 

,(?7cr) y Z^r] 



= {'VK"'', Vii''='^)^"^5™(if"'= - K'"^ - K'"= + K'"^), 

whence: 

-8i?3 ^ Eaii { (Vii^''^ VK'"^)iK'"' - K'"' - K^" + K'"^) 

■ {aauPcrabsPdvS'^^ -aauPcrPbsadvS"'' S''" - PauacrabsPdv^"'' S""" + PauacrPbsadvS"" S'"'')} 
= Ea6cd{(("a'"b)(^c,/3d) - {aa, Pb){Pc, ad) - {Pa, ab){ac, Pd) + {Pa, Pb){ac, ad)) 

■ {VK"\ \/K'^){K'''' - K"'^ - K'"' + K^^)} 

Relabeling the indices in the above expression yields: 

-8i?3 = J2abcd [^{aa,ab){Pc,Pd) - 2{aa,Pb){Pc,ad) - 2{Pa,ai,){ac,Pd) 

- 2{aa, Pb){Pd, ao) - 2{pa, ab){ad, Pc)] (V/^<^^ VK^'')K-'' 
= EabcdK'"'[Salab{VK->>)^WK'^''PJPc - 2alPb{VK-^)'^S/K^''alP, - 2plab{VK-^)^VK^''Pla, 

- 2alPb{VK'^YVK^''pJa, - 2p^ab{VK'^YVK'=''aJp,] 

= EacK'"'[SaliDa^°YiDp'^nPc - 2aliDp^'>YiDa^nPc - 2p^{Da^°Y{Dp^"^)ac 

- 2aUDP^r)^{Dp'^°'^)a, - 2p^ (Da^.n^ {Da^,°^)P4 

= Eac [^{Da'r ■ aa, Dp'r ■ Pc) - 2{Da'r ■ Pa + Dpl°^ ■ aa, Da'^'^' ■ p, + Dp^°'^ ■ a,)] . 
= Eac K''' [&{Fa{a, a),F,{p, /?)) - 8{Fa{a, P),F,{a, p))] , 

which is precisely (44). Alternatively, this can be derived from formula (40). 
• Computation of R4. It is the case that: 



20 



By Lemma 10: 

V (a B - B a ) n'^^^^^ ^'''~'> ^ (a B - B a )d K^°- (K^'' - K"-") 5^'^ 

Z_^aucr\'^auycr Hauler jy Z-^ aucr ^ yau^cr )^r^^ V 

= Ea„ {"a„ [ErdrKi^'iEcK^'M] - ^au [Er drK^^ {Z^K^^ Per)] 

- Pau [HAK^'iUcK^'o^cr)] + Pan [E.^.i^^"^ (Ec^"'"-)] } ^5^" 

= Eau {«a« [(VX«^ - (Vi^^^ /3h°0] - [(Vi^^^ - (Vi^«", a'r)] } <5^" 
So if we define the matrix Hia := Ylib \p°'^{P)(^bi ~ C"'^(a)/56i] , « = 1, . . . , -D, a = 1, . . . , we have: 

^/s ^r/ A/^ u=l C,7?=l 

Alternatively, this can be derived from formula (41). □ 

The denominator (34) of sectional curvature for (MP) is given by the simple formula: 
Proposition 11. For any pair of cotangent vectors a,/3 S T*C, 

abed 

Proof. Using double-index notation we may write equation (34) as follows: 

= Eabcd at, Per Pdv (K'^bgus j^cd^rv _ j^adguv j^bcgsr^ 

= Eabcdi^a^ "6) (/3c, Pd)K^'K^' ' Eabcdi^^a, Pd){ab, Pd)K-''K^^ 

and (51) follows by relabeling the indices. □ 
4.4 The rotationally invariant case 

Finally, suppose the Green's function K is rotationally invariant, i.e. that (K4) holds: 
K{x) = 7(||a;||), x , with 7 € C^{[{),oo)). 

We will use the convenient notation: 70 := 7(0), 7a6 := 7(||5° — 9^||), 7ab := 7'(||9° — ^''ll), and 
Tab '■~ i'^Wq^ ~ I^W) for a,b = 1, . . . ,N. Then we can evaluate the first and second derivatives of K: 

Lemma 12. For rotationally invariant kernels, it is the case that: 

VK(x)=7'(l|a^ll)^, (52) 
r,2Trf^\ L"nuih 7'(l|a;||)l a;a;^ ^ 7'(l|2;||) ^ 

D K{x) = [7 (Ikll) - + -j^ In (53) 

-7 UI^IIJ||^||2 + ii^ii W' 

where Id is the DxD identity matrix and Pr"'"(x) := Id — is projection to the hyperplane o/R-^ 
that is normal to x. 
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Proof. We have that diK{x) = 7'(||a^||) and (52) follows immediately. Also, 

5.Wa;) = MsfT7'(l|x||)+7'(ll^ll)^z|7:^^+7'(ll^liyz|T^ 

= 7"(I|:^II)S; + ^.5- - YdklDgp = [Y'dklD - ^] gp + ^<5-, 

which implies (53). □ 

Because of (52), in the rotationally invariant case, the "scalar compression" C^ia) really does 
measure a multiple compression of the flow a" between q°- and q^. We can decompose the vector 
strain S"'''{a) into the part parallel to the vector q°- — and the part perpendicular to this: let 

S'"*'(a)ll := {S"-\a), u"^), S''\a)^ := S''\a) - S''\a)^\ u^K (54) 

Note that S'°^(q;)" is a scalar while S'^^{a)-^ is a vector. In particular it is the case that C°^(q:) = 
7^^-S'"*'(a)l'. Moreover, formula (53) allows us to simplify the first term Ri in the curvature formula. 
Substituting (53) into (42), we get the rotationally invariant case for Ri: 

Proposition 13. In the rotationally invariant case (K4) we have that 

Ri =E (^(5"'(a)ll Pa - S'^W aa, 5«^(a)ll pt, - 5«''(^)ll a,) (55) 

ajtb 

^ab / Qdb / ^ Q co-b/p^J- ^ ^, Qabt^\l. ^ p oab ( q\1- 



+ Oil a ' 611 (^°''(^)^ ®Pa- S'^\ii)^ ® aa, S-\a)^ ® /Jfe - S'^\p)^ ® a,)) . 

9 II ' 
In the above we use the inner product of tensor products, {vi ®wi,v2® W2) ■= {vi,V2){wi,W2). 

Proof. For any pair of covectors t] and /z in T*£, by (53) we have that: 

S''\rif D'^K"-^ S''\n) = i^i,S"-\f^fv"'' {u"''f S'^^fi) + ^^^j'^^^ S^^-qf Pr^ {u"'') S"\n) 
= l^bS'^'ivy^ S-\n)\\ + (5»''(7?)^,5«''(m)^). 
Inserting this expressions into (42) yields the desired result. □ 
4.5 One landmcirk with nonzero momenta 

A simple special case is when only one landmark carries momentum. We now compute the numerator 
of sectional curvature when both cotangent vectors are nonzero at only one of the £)- dimensional 
landmarks {q^ , . . . ,q^). We define: 

{T;C)i := {tj e T*C I rya = for a > 1} 

so that the elements of the above set are cotangent vectors of the type = (771, 0, . . . , 0). 
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Figure 3: Dragging effect of one momentum-carrying landmark (bullet •) on a grid of landmarks 
(circles o), with 7(x) = exp(— i^), a = 1.5. Left: initial configuration, with initial momentum 
Pi — (2.7,1.8) also shown. Right: configuration after one unit of time, with trajectory of also 
shown; the grid represents the diffeomorphism <fiQi, obtained by integrating a'^" in time. 



Proposition 14. In C^{R^), for any pair a, /S e {T*C)i the four terms o/ i?(a'*, /S", /?*•, a") are 
given by Ri ^ R2 ^ R3 ^ and R^ = -| E^b=2(^''' ■^^^)r" (K~^)a6, where 

H"" :- (7ai -7o)((ai, VX"i)/3i - VX'^i)ai) , for a > 1. 

Proof. The vanishing of Ri can be checked directly (note that the sum in (42) is taken over a ^ b 
since S'^''{ri) = for all c and r/). Also, using formula (39) we see that all mixed forces Fa are zero, 
therefore i?2 = i?3 = by formulae (43) and (44). Also, by (41), D°(a, (3) = (7ai -7o)(ai, VK°i)/3i 
since a, /3 G {T*£)i; a similar expression holds for D°-{f3, a), which concludes the proof by (45). □ 

Therefore when a,/3 Cz {T*C)i the sectional curvature is always negative; we can understand 
this by considering the geodesic flow in this case. It follows immediately from Proposition 6 that if 
we start with zero momenta Pa at all a > 1, then the momenta at these points stay zero, while 
the momentum at q^ remains constant. Thus the velocity of q^ is just given by K{{))pi and this is 
constant. The point q^ carrying the momentum moves in a straight line at constant speed, while 
the other points 5° (a > 1) are carried along by the global flow that the motion of q^ causes and 
move at speeds = /sT'^^pi, which are parallel to q^ (but not constant). As shown in Figure 3 (the 
central landmark q^ is the only one carrying momentum) what happens is that all other landmark 
points are dragged along by q^, more strongly when close, less when far away. Points directly in 
front of the path of q^ pile up and points behind space out. 

Negative curvature can be seen by the divergence of geodesies. If you imagine slightly changing 
the direction of pi in Figure 3, the final configuration of the landmark points (say, after one unit 
of time) will differ greatly from the one caused by the original value of pi. Also, if you imagine q^ 
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moving along two nearby parallel straight lines, the differential effect on the cloud of other points 
accumulates so that the final configurations will differ everywhere; thus, even though the initial 
landmark configurations are close, the final configurations will be far away. In general, the last 
negative term in the curvature expresses the same effect: the global drag effect of each point results 
in a kind of turbulent mixing of all the other points (think of a kitchen mixer the motion of whose 
blades mixes the whole bowl). 

Proposition 14 simplifies in the case of £ = £^(M^) (two landmarks only). We shall write: 

4~{a^y^), ai:=a^-4u^\ p\ := {p^u^^) , fit := - (56) 

Proposition 15. In the case of C = £^(K-°), when a,(3 G {T*£)i the numerator and denominator 
of sectional curvature are given by, respectively: 

R{a\p\p\a^) = i?4 = 70 {l[,f \\p\ai - af/3^f , (57) 

4 7o+7i2 

Mt^c WWl'C - P?T'C = 70 - + (58) 

Again we have used the inner product of tensor products {vi ®w\,V2® W2) ■= {vi,V2){wi,W2)- 

Proof. It is the case that R4 = — |||i/^|p(K^^)22 (matrices H"^ were defined in Proposition 14). 
But (K~^)22 = (7o — 7i2)~^70i whereas from Proposition 14 we have 

WHY = (70 - 7i2)'((ai, V/^^2)2||A||2 + (A, V0'll«if - 2(ai, VO(A, VO(«i, A)), 

where VK^"^ = 7^2^^^ by (52). Inserting expressions (56) into the above formula yields (57). 
From Proposition 11 we have that the denominator is given by 7o(||Q!i|P||/3i|P — (ai,^i)^); again, 
inserting (56) into such formula yields (58). □ 

We will generalize the above results in the next section. 



5 Landmark geometry with two nonzero momenta 

The complexity of the formula for curvature reflects a real complexity in the geometry of the land- 
mark space. But there is one case in which the geometry such space can be analyzed quite completely. 
This is when there are only two nonzero momenta along a geodesic. To put this in context, we first 
introduce a basic structural relation between landmark spaces. 

5.1 Submersions between landmark spaces 

Instead of labeling the landmarks as 1. 2, • • • , N, one can use any finite index set A and label the 
landmarks as q"' with a € A. And instead of calling the landmark space £^ , we can call it 
Now suppose we have a subset B C A. Then there is a natural projection tt : jC-^ CP gotten by 
forgetting about the points with labels in ^ — B. In the metrics we have been discussing this is a 
submersion. In fact, the kernel of d-K, the vertical subspace of T£-^, is the space of vectors v"' such 
that v°- = it a G B. Its perpendicular in T* is: 

{T*C-^)b := {p e r*£-^ \pa=OioJ: a&A- B) 
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so the orthogonal complement of ker(d7r) in TC/" is the space of vectors where p is in (T*£'^)b. 
On this subspace, the norm is just 

6,6' eB 

whether is taken to be a tangent vector to A or to B. In other words, the horizontal subspace 
for the submersion tt is the subbundle (r*£-^)g C T£-^ of tangent vectors where p has zero 
components m. A — B and this has the same metric as the tangent space to CP . In particular, 
from the general theory of submersions, we know that every geodesic in CP beginning at some 
point 7r({(7''}) has a imiquc lift to a horizontal geodesic in starting at {<?"}. The picture to have 
is that all the landmark spaces form a sort of inverse system of spaces whose inverse limit is the 
group of diffeomorpisms of MP . 

We don't want to pursue this is in general, but rather we will study the special case where the 
cardinality of B is two. We might as well, then, go back to our former terminology and consider the 
map TT : — )■ C? gotten by mapping an iV-tuple {q^,q^, • • • , g^) to the pair {q^,q^). Moreover, we 
want to consider only the case in which the kernel K is rotationally invariant as in (K4). A basic 
quantity in all that follows is the distance p := \\q^ — between the two momentum bearing points. 



5.2 Two momentum geodesies 



Remarkably, we can describe, more or less explicitly, all the geodesies which arise as horizontal 
lifts from this map. These are the geodesies with nonzero momenta only at and q^. Moreover, 
the formula for sectional curvature for the 2-plane spanned by any two horizontal vectors can be 
analyzed. This analysis was started in the PhD thesis of the first author [23] and has been pursued 
further in [22]. 

The metric tensor oi jC = in coordinates is obtained by inverting the 2x2 matrix K: 



K = 



70 


lip) 




. lip) 


70 





(K-i)i2 



(K-l)22 

(K-i)2i 



il' 



o-lip)')-'lo 

iii-iip)')-'iip) 



(59) 



so that the cometric and metric, for all covectors a, /3 £ T*C and vectors v,w G TqC, are simply: 



5-i(a,^) =7o((ai,/3i) + (a2,^2)) +7(p)(("i,/02) + (a2,/Si)), 



(60) 



g{v,w) = 



ll 



lip)^ 



The geometry of the two-point space is best understood by changing variables for the landmark 
coordinates (g^,g^) and the momentum (pi,P2) to their means and semi-differences, that is: 



q ■■= 

,1 



2 ' 

so that: q'^ =q + Sq, 
Then the cometric (60) becomes: 



6q := 

„2 



q 



2 

q-6q, 



p := 

Pi 



Pi +P2 

2 ' 
p + dp, 



5p := 

P2 = 



Pi -Pl 

2 ' 
p — Sp. 



g-' ((«, da), (p, 513)) = 2(70 + 7(p)) («, P) + 2(70 - lip)) {5a, 



(61) 



With these coordinates, the two-point landmark space becomes a product ^ x in which all fibers 
V X {5qo} are flat Euclidean spaces though with variable scales, all fibers {gg} x Vs are conformally 
flat metrics sitting on the manifold — {0} and the tangent spaces of the two factors are orthogonal. 
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Proposition 16. In terms of means and semi-differences, the geodesic equations for C"^ {MP) are: 
%={lo+ lip)) P, P = 0, 

k = (70 - 7(P)) Sp = -2^ (llpf - \\Spf) 5q. ^^^^ 

The above result is proven by direct computation. We can solve these equations in four steps. 

1. First the linear momentum p is a constant, so "center of mass" q moves in a straight line 
parallel to this constant: 

q{t) = m + ( ^ (70 + l{p{r)))dT) p. (63) 

2. Secondly, if we treat vectors 6q and 6p as 1-forms in M.^, equations (62) also show that: 

(Sq A dp)' = Sq A 6p + 5q A 5p = [(scalar) 6p] A6p + Sq A [(scalar) 5q] = 0, 

so the angular momentum 2-form 6q A6p G /\^ is constant; we write this as w A where co 
is the nonnegative real magnitude of the angular momentum and (e^,e^) is an orthonormal pair. 
Then it follows that: 

Sq{t) = ^p{t)[cos {9{t))e^ + sin {9{t))e'^], for some function e{t). 

3. Thirdly, we can express 9{t) as an integral: 

6q= 5p[cos(6')ei +sin(6l)e2] + ^pe[ - sm{e)e^ + cos{0)e'^], so 
5q ASq = -\p^0e^ A e^ as well as (from (62)): 
SqASq= (70 - l{p))5p A Sq ^ -uj (70 - 7(p)) A e^; 
combining the second and third lines, wc find: 

eit) = ei0)+4.f ^°-j\f^^ dr; (64) 
Jo P{t)^ 

note that by (Kl) and (K2) it is the case that 70 > j{p) for any p > 0, so ^ is a monotone increasing 
function if w 7^ 0, otherwise it is a constant. 

4. The last step is to solve for p{t). This can be done using conservation of energy [16, p. 51]. 
Equations (62) are in fact the cogeodesic equations for the Hamiltonian 'H(p, q) of section 4.1, which 
we may rewrite in terms of means and semi- differences as 

H = (70 + 7(P)) llpf + (70 - 7(P)) 
by (61); hence this function of p and \\6p\\ is a constant [p is also a constant). Then we calculate: 

ip^y = 4{Sq,Sqy = 8{Sq,Sq) = 8(70 - 7{p)){Sp,Sq) =^ p = A^^^^M {6p, Sq) . 



But: 



V2 , ,.2 /r„ .„X2 , ,|c . r 1,2 „ri|2nr„n2 P^ f ^ - {-fo + l{p))\\p\\' 



{Sp, Sq)' +uj' = {Sp, Sq)' + \\Sp A SqW = ■ IM' = 



4 V 70 - 7(P) 

I - ^p2 [n ^ (^0 + lip)) IIpIP] - 4a;^ (70 - lip))- 
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Figure 4: Converging and diverging trajectories for two landmarks in two dimensions. In these ex- 
amples ^(x) = exp(-ix2), (gi(0), g2(oj) = ((1, 0), (-1, 0)), (pi(0),p2(0)) = ((-10, 8.6), (10, -8.6)) 
for the graph on the left, (331(6), 332(0)) = ((-10,9), (10, -9)) for the graph on the right. The thick 
and thin lines are, respectively, the trajectories of the first and second landmarks. 

This means that the function p(f) is the solution of: 

t = / -y=, where: F{x) := Hx'i'yo - 7(x)) - Wpfx^d'o ' li^f) ' ^^'ilo - 7{x)f. 
Jp{o) 2^F{x) 

(65) 

Summary. If we fix constants H, p, w, p{0), 0{O), q°'{0) (for all a), we can first integrate (65) 
to get p{t) (the separation of and q"^), then integrate (64) to find their relative angle 9{t), then 
integrate (63) to get their center of mass q{t). This gives the trajectories of q^ and q^. The remaining 
points are dragged along as solutions of: f^q^it) = 7(||g''(t) - q^{t)\\)pi{t) +7{\\q"{t) - q'^ {t)\\) P2{t) . 

As worked out in [22], one can classify the global behavior of these geodesies into two types. One 
is the scattering type in which q^,q'^ diverge from each other as time goes to either ±00. This occurs 
if the linear or angular momentum is large enough compared to the energy. In the other case where 
the energy is large enough compared to both momenta, they come together asymptotically at either 
t — +00 or —00, diverging at the other limit. In both cases, they may spiral around each other an 
arbitrarily large number of times (see Figure 4). 

5.3 Decomposing curvature 

Next we consider £^(M^): we want to compute the sectional curvature R{a\f3^, /?", a") for cotangent 
vectors that are nonzero at only {q^,q'^). Also, we will use the notation u := ||^i~^2|| for the unit 
vector from q'^ to q^ as well as p = — for their distance. Similarly to (54), we will also want 
to decompose any vector in r/ € M.^ into its parts tangent to u and perpendicular to u: 

ry" := (ry, u), and r?""" := 77 — r?" u. 
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Once again note that 77" is a scalar whereas 77^ is a vector. Following the notation used to describe 
geodesies above, for any a S {T*C)x^2 '■= {^7 G T*C | rya = for a > 2}, we write a = \{a\ + a^) 
and 5a = 5(0:1 — 02) • 

Proposition 17. In £^(K^) for any pair a, jS € {T*C)i^2, the terms i?i,i?2 and R3 in the numer- 
ator of sectional curvature can be written as 

Ri = 4(70 - 7(P)) V(P) - <5/3'lai, fo'l/?2 - Sl3\\a2) 

+ 4(70 - 7(P))^ {Sa^ ®fii- 5fi^ (g> ai,Sa^ (E> p2 - Sp-^ (g) 02), 

i?2 = -4(70 - 7(p)) 7'(p)2 (^all/3i - ^^llai, (5q;II/?2 - ^/3IIq2), 
i?3 = ^" y(p)^[((ai,/32) + (A,a2))' -4(ai,a2)(A,fe)]. 

We need the following result. 

Lemma 18. For any a € (T*£)i_2, the discrete strain S^^{a) is given by: 

S'\a) = 2{-/o-^{p))Sa. (66) 
For any pair a,j3 & {T*C)\^2 it is the case that Fa{a, /3) = for a > 2, whereas 

F,{a,p) = -F2{a,p) = ^^{{ax,^^) + {Pi,a2))u. (67) 

Also, 

D\a, p) = 2(70 - 7(P)) i{p) 5a\\ p2, D^ia,p) = 2(70 - j{p)) i{p) fo" ft. (68) 

Remetrk. We are not interested in D'^{a, (3) for a > 2 since the terms in formula (43) where they 
appear are zero (because Fa{a, /?) = for a > 2). 

Proof of Lemma 18. The formula for the discrete strain results from: 

5^2(a) = {a^f - {Jf = E6(-fs^" - ^^'')a6 = 70^1 + 7(^)^2 - 7(p)ai - 70^2 = 2(70 - l{p))5a. 
The values for F follow immediately from formula (39) and WK^^ = j'{p) u. Note that: 

C'^a) = C^\a) = {S''{a),VK^^) = (2(70 - j{p))San\p)u) = 2(70 - 7(p))7'(p) so 

£)i(a,/3) = Ci2(a)^2 = 2(7o-7(p))7'(p)^a"^2 and Z?2(c,,^)=c2i(a)/3i=2(7o-7(p))7'(p)'^a"/3i. □ 

Proof of Proposition 17. The Ri expression follows by substituting the expressions in (66) into for- 
mula (55), noting that the only non-zero terms in the latter are for (a, b) = (1, 2) and {a, b) = (2, 1). 
By Theorem 9 and the fact that F2 = —Fi from Lemma 18, R2 is given by 

R2 = {D\a, a) - D\a, a),Fx{P, /?)) + {D^ (/3, /3) - D^p, /3), Fi(a, a)) 

- {D\aJ) - D\a,(3)+D\(3,a) - D\(3,a), Fx{a, (3)). (69) 
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Again by Lemma 18 wc have that D^{riX) — -D^(?7, C) = ^4(70 — jip))^' (p) Stj^^ for any pair 
r?,C G (T^*£)i,2, while Fi{r],C) = |7'(p)((»?i, C2) + {m,Ci))u- Applying this to all the terms we get 
the expression for R2 in the statement of the proposition. 
As far as Rs is concerned, by Theorem 9: 

= 70 /3),Fi(a, /3))- a),Fi(/3, /3))] +7(p) [(Fi(a, /3),F2(a, /3)) - (Fi(a, ^ 

+ 7(p) [{F2ia, f3),Fiia, /?)) - {F^ia, 13))] +70 [{F^ia, P),F2{a, - (i^2(a, a),F2{P, (3))] 

= 2(70 - 7(/3)) [{Fi{a, F^ {a, /?)) - {F, (a, a),F, /?))] 

= ^^i^(7'(p))'[((ai,/32) + (/3i,a2))^ -4(ai,a2)(/3i,/32)], 
where we have used the fact that F2 — —Fi, by equation (67). This completes the proof. □ 

The expressions provided by Proposition 17 become much clearer if we go over to means and 
semi-differences, i.e. if we use the substitutions: 

ai=a + Sa, a2=a-Sa, Pi=P + 6P, P2=P-Sp. (70) 

Corollary 19. For any a,l3e (T*£)i,2, with C = £^(]R^), it is the case that: 

i?i = 4(70 - 7(/5))%"(p) (||<5^lla - SJM^ - \\Sl3ha - Sahl3\\^) 

+ 4(70 - 7(p))'^(||<5/3^ (g,a-Sa^(g, - \\Sl3^ (g, Sa - 6a^ (g Spf), 

R2 = -4(70 - 7(P))7'(P)' (||(5/3"a - ,5all^||2 - \\5pha - dahpf), 

R3 = (70 - 7(p))7'(p)^ {2\\S/3(g)a-Sa(g)]3f - \\]3(g)a-a^pf-\\S/3^da-Sa®6l3f). 
Proof. By insertion of formulae (70) it is easily seen that 

{5Jf3i - 5/3" ai, foil /32 - (5/3lla2) = \\S(3^^a-Sa^f - \\5pha - 6Jd(3\\^ , 

so the new expressions for Ri and R2 follow immediately. Also: 

[(ai,/32) + (/3i,a2)]^ -4(ai,a2)(/3i,/32) 

= [2{{a, 'P) - {5a, 5(3))] ^ - 4((a, a) - {Sa, Sa)) {(fi, 'p) - {5(3, 5(3)) 

= -2[2{{a,a){~p,'p) - {a,W)\ - 2[2{{5a,5a){5f3,5f3) - {5a,5p)^)] 

+ A[{a,a){5f3,5f3) + {P,T3){5a,Sa)~-2 {a, T3){5a,5l3)] 

= -2\\Ji®a-a®]if -2\\5l3®5a-5a®5l3f +A\\5l3®a-5a®^\\^. □ 

The fourth term i?4 is the only one which involves the other points q°',a > 2. But one has an 
inequality for this term involving the same expressions in a and /3: 

Proposition 20. Any pair a, (3 E (T* C^)i^2 are constant 1-forms on which are pull-hacks via 
the submersion — > of constant 1-forms on . We can therefore consider the curvature term 
Ri{C^) = -|||M,/3"]£"f on and the corresponding term Ri{C?) = -|||M,^'']£2||2 on C? . 
Then we have the inequality: 

\ 2 



R^C"") < Ri^C) = -67' (p)' 



^^!^-^^^^||5^lla-foll/3||2 + (7o-7(p)) \\6pha-6a\ 
70 + 7(P) 
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Figure 5: Typical covectors a = (ai,a2) in spaces JUt*/;, (J-^T*/;, and T*/;, for £ = /;2(k2)_ 

Proof. Firstly, note that [a'j^'J^w breaks into perpendicular parts: a vertical part in the kernel of 

diT and a horizontal part which is simply the horizontal lift of [a\ (i'^]c'^ ■ This explains the inequality 
assertion in Proposition 20. To calculate i?4(>C^), we use the last expression in (45), i.e. 



70 + 7lPj L 

- 27(p)(fcll/3i - ^^llai,(5all/32 - ^/3"a2)| 



where we have used (59) and (68). The final result follows after inserting (70) into the above 
expression and performing some algebra. □ 

Note that all terms in Corollary 19 and Proposition 20 are very similar. In fact, they are all 
"components" of the norm Ha A of the 2-form whose sectional curvature is being computed. 
First note that we can decompose T* into the direct sum of three pieces, namely: 



|(au, — a?i) I a G M}, 

{(p, -p)|peM^,p±M}, 
{(p,p)|peM^}, 



dim((5llT*£2) = 1^ 
dim(5-Lr*£2) =£)-!, 
dim(r*£2) 



where as usual u := n^i^^ay (sec Figure 5). Note that these throe subspaces are orthogonal with 
respect to the cometric by virtue of (60). An arbitrary covector a = (ai, a2) G T*C^ can be uniquely 
decomposed into the summation a = a(i) + a(2) + Q^(3)) with: 

«(!) := {SJu, -SJu) G S\^T*C^, a(2) := (fo-^, -Sa-^) G 5^T*C'^, a^) := (a, a) e T* C"^ . (71) 

So it is the case that: (i) a € (5llr*£2 ^ J-Lq, = and 5 = 0; (ii) a G J-Lr*/:^ ^ = and a = 0; 
(iii) a e T*£2 ^q,|| = q and ^a-^ = 0. 

Consequently the space of 2-forms /\^ T*£^ decomposes into the direct sum of five pieces: 



/\r*£2 = 0y,, with: 

i=l 

^2 := ^^T*£2 A T*£ 

2 

Vi := /\{S^T*£^), 



Vi := (5"T*£2 a r*£2, 
:= 5llT*£2 A ^-LT*£2, 

^5 := Ar-c')- 
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(Since (5llT*£2 

is one-dimensional it creates no 2-forms.) Once again, note that the spaces Vi, . . . , V5 
are pairwise orthogonal with respect to the inner product 

(q A ^, ^ A ri)^2 := {a, Or^r^ ('^' ^)t*£2 " ^)t'C^ Ot-c^ ' /5> C> e T*C'' (72) 

by the orthogonality of 6^^T*£., 5^T*C, and T*£. Any 2-form a A /3 then decomposes into the sum 
of its five projections onto these subspaces and its norm squared is the sum of the norm squared of 
these components. Let us first give the five pieces of its norm names: 

Ti := ||(5^"u(^a- (5a"u(g)/3||^ 
T2 := \\6I3^ ®a-5a^ T3 := ||(5/3llu fo-^ - foUu (g) J/J-^f, 

Ti := \\5^^ ®5a^ -5a^ ®5p^\\^, T5 := ||^ a - a (g) ^||^ 

In the above definitions || j| indicates the Euclidean norm. We have to be careful here: we have 
been using Eiiclidcan norms in in all our formulas above and now we arc dealing with norms in 
T* C^-. these essentially differ only by a factor, by (61). More precisely, the following result holds: 

Proposition 21. The denominator of the sectional curvature (17) /or £^(R^) can he written as: 

\\a A = 4(72 - 7(p)')(Ti + T2) + 2(70 - 7(P))'(2T3 + T4) + 2(70 + ^{p))^n. (73) 

Proof. We may apply decomposition (71) to both a = X]f=i (^{i) ^^'^ P = /^(j)) and write 

aAP =(«(!) A /3(3) - A a(3)) + (a(2) A /3(3) - ^(2) A a(3)) 
+ ("(1) A /3(2) - A q;(2)) + a(2) A /3(2) + a(3) A ^(3), 

where the five summands on the right-hand side belong to Vi, . . . , V5 respectively. We have 

A /3(3) - A a(3)||^2y.£2 = 
= ||a(i) A /3(3)ll^2 j,.£2 + ||/3(i) A Q;(3)||^2 j,.£2 - 2(a(i) A /3(3),^(i) A a(3))A='T*£2 

= 4(7o' - lipf) [{Sa^mf + {SP^^rm' -26JSl3Ha,P)] = 4(70^ - ^ipf) T,, 

where we have used (72) and (61) in step (*). The square norm of the remaining four terms is 
computed similarly. Orthogonality of Fi, . . . , V5 finally yields (73). □ 

To express the formulas for the numerator of sectional curvature succinctly, let us also introduce 
abbreviations for the coefiicients involving 7: 

fci(p):=(7o-7(p))V(p), Hp) 

k3{p)--={lo-lip)h'{pf, h{p) 

Note that ki, fc2, ^3 and k4 are all homogeneous of degree 3 in 7 and degree —2 in the distance p 
or dp on £^ . Moreover ^2 is negative, ^3 and ki are positive, while fci may be positive or negative. 
For all 7 of interest, 7' is everywhere negative, starting at decreasing to a minimum at some po, 
then increasing back to at 00. Then ki is negative for p < po and positive for p> po. 
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(70 - 7(P)) 
(70 -7(P))' 



27 [p) 



(74) 



70 + 7(P) 



I'ipf- 



Figure 6: The coefRcients of Ti (top left), T2 (top right), T3 (bottom left) and T4 (bottom right) for 
the Bessel kernels 7 (shown with thin lines) and the Gaussian kernel (shown with the thick line). 
The kernels are scaled to normalize 7(0) and 7"(0). 

The following equalities are proven by direct computation: 

\\Sp^a-Sa^W = Ti+T2, 
\\6l3-^ (g)Sa- 5a^ (g) Sl3\\^ + T4, 
\\SP (S)6a-Sa® Sf3\\^ = 2T3 + T4. 

Inserting notation (74) and the above equalities into Propositions 17 and 20 immediately yields: 
Proposition 22. We can write the terms in the numerator of sectional curvature for C^{R^) as: 

i?l =4fci(Ti-T3) + 4fc2(T2-r3-T4), i?2 = -4fc3(Ti-T3), 

i?3 =A;3(2(Ti+T2)-2r3-T4-T5), = -GiksTs + kiTi), ^ ^ 

hence R = R{aK , a«) = J2Ui ^ 

may he expressed as: 

R = 2(2fci - fca - 3fc4)Ti + 2(2^2 + fc3)T2 + 4( - fci - fca - fc3)T3 + ( - 4^2 - ^3)74 - fcsTs. (76) 

By virtue of Proposition 20 the above proposition still holds in the case of £^{M.^) as long 
as a, /3 e (T*£)i_2 and the equality signs for R4 in (75) and R in (76) are substituted by "<". The 
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GEODESIC 1 

GEODESIC 2 

-0 2 I ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' 

0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 2 4 6 8 10 12 

P q' 

Figure 7: Left: sectional curvature JC for £^(M'^) (from Proposition 23), as a function oi p = Iq^ — q^l; 
here 7(0;) = exp(— ^x^). Right: two trajectories in £^(M^) shown in the {q^,q^) plane (under the 
assumption that < q^). Both geodesies originate at [q^^q^] — (0,4), and lie in the region 
where /C > (above the upper dotted line, that indicates the zero of K. at \q^ — q^| ~ 1.53); they 
have different initial momenta {pi,P2) = (1, 1) and {pi,P2) — (1,0.4), and exhibit conjugate points. 

coefficients in (76) may have all sorts of signs for peculiar kernels. However, the kernels 7 of interest 
are the Bessel kernels (3) and the Gaussian kernel, which is their asymptotic limit as their order 
goes to infinity. The coefficients for these kernels are shown in Figure 6. We see that the coefficients 
of T2 and T3 are negative while those of T-i are positive. Henceforth, we assume we have a kernel 
for which this is true. 

5.4 Sectional curvature of £^(]R^) 

Finally, we will now explore the important example of two landmarks on the real line. In this partic- 
ular case the manifold is two dimensional, so sectional curvature /C will turn out to be independent 
of cotangent vectors a and /3. In fact, given the translation invariance of the metric tensor, it will 
only depend on the distance p = \q^ ^ q^\ between the two landmarks. 

The spaces S'^^T* C'^{M}) and T* C'^{M}) are one-dimensional while 5^T*C:^{M}) = {0}. Thus 

2 

f\T*C^{M}) = 5llT*£2(Mi) A T*C^{M}) 

and the only non-zero term in (76) is Ti. Therefore combining formulas (73) and (76) we get: 
Proposition 23. The sectional curvature of C'^iJS}) is given by 

^ 2ki - fcs - 3fc4 ^ 7o - lip) II _ 27o~7(p) , .2 
2(7o' - lip?) 70 + lip) (70 + lip))' ■ 
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GEODESIC 1 

GEODESIC 2 

_gl 1 1 1 , 1 1 

-6 -4 -2 2 4 6 

Figure 8: Existence of conjugate points in £^(K^), with j{x) = exp(— Both geodesies 
originate at landmark set {q^,q^) = ((— 2, — 4), (2, 4)) ; the first one (dashed) has initial momen- 
tum {pi,P2) = ((0, fO), (0, 10)) e r*£^ while the second one (continuous) has initial momentum 
{pi,P2) = ((6, 10), (-6, 10)) e 5llT*£2 e T*£^. The geodesic trajectories exhibit conjugate points. 

The above function /C is shown on the left-hand side of Figure; 7 as a function of p, for the 
Gaussian kernel. The coefficient of the term Ti in (76) is negative for p small and positive for p 
large. The "cause" of the positive curvature has been analyzed in [23]. Roughly speaking, suppose 
two points both want to move a fixed distance to the right. Then if they are far enough away, they 
can just move more or less independently (we shall refer to this as Geodesic 1). Or (i) the one in 
back can speed up while the one in front slows down, then (ii) when the pair are close, they move in 
tandem using less energy because they are close and finally (iii) the back one slows down, the front 
one speeds up when they near their destinations (Geodesic 2). This gives explicit conjugate points 
(in the sense that two points arc joined by distinct geodesies) and is illustrated on the right-hand 
side of figure Figure 7 (where Geodesies 1 and 2 are represented, respectively, by the dashed and 
thick curves). 

5.5 Sources of positive curvature; obstacle avoidance 

There is another source of positive curvature in £^ in higher dimensions. It is clear from equation (76) 
and Figure 6 that any positive curvature must come from the term with Ti or the term with T4. As 
the five terms are orthogonal, we can make all of them but one zero. 

For example, if we choose a = {Sa^^u, —Sa^^u) G (5llT*£ and /3 = (/3,^) e T £, then it is the case 

that Ti = ((5all)^||/3|p and it is the only non-zero term. Then, if p is sufficiently large, the sectional 
curvature for this 2-plane is positive as discussed in the last section. Figure 8 illustrates an instance 
of the existence of conjugate points for two geodesies in £^(R^); the momenta {pi,P2) of each of the 
two trajectories belong at all times to (JUtX^ ® T*£2. 
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The other possibihty is that T4 is the non-zero term, which happens when a = ((5q;^, —Sa-^) £ 
5^T*C and ^ = {513^,-513^) e 6^T*C. We have = 2{\\6a^f\\5l3^\\'^ - {5a^,5l3^f), and for it 
to be nonzero it is required that D >2> because T4 is the norm of a 2-form in A {5^T*C), which 
has dimension (D — 1) (D — 2)/2. The positive curvature of this section is readily seen by considering 
the geodesies which these vectors generate. The simplest example is the following: 

Proposition 24. The circular periodic orbit of radius r: 

q\t) = {rcost,rsmt), q^(t) = -q^{t), (77) 

t gB., is a geodesic m£^(]R^) if and only ifr is the solution of the equation ^o — ^{2r) {2r) = 0. 

Proof. For orbit (77) it is the case that q = 5q = q^ and p = 2r; also pi ~ (70 — 7(p)) ^q^ and 
P2 = —pi , so that p = and Sp = pi. The first three equations of (62) can easily be checked; the 
fourth one holds if and only if 70 — 7(2r) + r^'{2r) =0. □ 

(The above result was also proven by Frangois-Xavier Vialard of Imperial College, London.) 
Orbit (77) has the property that at time tt, q^ and q^ interchange their positions: it is a geodesic 
from the set of landmark points ((r, 0), (-r, 0)) e ^^{M?) to the set ((-r, 0), (r, 0)) G C'^{M?). But 
if these points live in K^, they can move around each other in any plane containing the points. Thus 
we have a circle of geodesies in £^(K^): 

q^{t) = {r cost,r cosO smt,r saiO smt), ffit) = —<l^{t) 

all connecting ((r, 0, 0), (-r, 0, 0)) to ((-r, 0, 0), (r, 0, 0)) , for any 6 e [0, 27r). This is exactly like all 

the lines of fixed longitude connecting the north and south pole on the 2-sphere and means that one 
set of landmark points is a conjugate point of the other in £^(K^). This is the simplest example of 
how geodesies between landmark points must avoid collisions and so make a choice between different 
possible detours, leading to conjugate points and thus positive curvature. 

6 Conclusions 

We believe that £^{R^), the Riemannian manifold of A'' landmark points in D dimensions, is a 
fundamental object for differential geometry and that we have only scratched the surface in its study. 
We started with a basic formula which computes sectional curvature of a Riemannian manifold in 
terms of the cometric, its partial derivatives, and the metric itself (but not its derivatives). This is 
particularly adapted to computing curvature for manifolds which arise as submersive quotients of 
other manifolds and gives O'Neill's formula as a corollary. We then applied this to derive a formula 
for sectional curvature of the space of landmarks. This formula is not simple but, like Arnold's 
formula for curvature of Lie groups under left- (or right) -invariant metrics, splits into a sum of four 
terms. The four terms involve interesting intermediate expressions in the two vectors (or co-vectors) 
which define the section and which have relatively simple geometric interpretations. We called these 
the mixed force, the discrete vector strain, the scalar com,pression and the landm,ark derivative. The 
geodesic equation in its Hamiltonian form is quite simple and involves the force as expected. We 
also gave several concrete examples to illustrate the nature of these geodesies. 

Finally, we have examined in detail the case of geodesies in which only one or two landmark points 
have non-zero momenta, and computed the curvature in sections spanned by such geodesies. We 
found that in this case there are essentially two sources of positive curvature. One can understand 
them through the non-uniqueness of geodesies joining two A''-tuples: the first sort of non-uniqueness 
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is caused by the two points with non-zero momentum choosing between converging in the middle of 
the geodesic ( "car-poohng" ) or moving independently and not converging; the second occurs only 
when D > 3 and arises when the same two points need to get around each other and must choose 
on which side to pass (if D = 2, this sort non-uniqueness also occurs but comes from non-trivial 
topology, not curvature). 

One of the most important questions left open is explore how prevalent positive curvature is in 
general, i.e. for geodesies in which all points carry momentum. Answering this question is central 
to applications of landmark space in which geodesies are actually computed. One might hope that 
the picture for two momenta is true in general but this is far from clear. It seems interesting to 
explore whether there is some sort of "index" for curvature forms — a numerical measure of how 
much positive vs. negative curvature is present. Another important question is to explore the shape 
of the coefficients in (76) for different kernels. More generally, what is the impact of different kernel 
types (Bessel, Gaussian, Cauchy) on the corresponding geodesies? Finally, note that all kernels have 
a length constant built into their definition so the geometry of the space of landmarks is far from 
scale invariant. Thus one should analyze what happens asymptotically when the points are very 
close relative to this constant or are very far from each other. 
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This paper deals with the computation of sectional curvature for the manifolds of A'^ land- 
marks (or feature points) in D dimensions, endowed with the Riemannian metric induced by the 
group action of diffeomorphisms. The inverse of the metric tensor for these manifolds (i.e. the 
cometric), when written in coordinates, is such that each of its elements depends on at most 2D 
of the ND coordinates. This makes the matrices of partial derivatives of the cometric very 
sparse in nature, thus suggesting solving the highly non-trivial problem of developing a for- 
mula that expresses sectional curvature in terms of the cometric and its first and second partial 
derivatives (we call this Mario's formula). We apply such formula to the manifolds of landmarks 
and in particular we fully explore the case of geodesies on which only two points have non-zero 
momenta and compute the sectional curvatures of 2-planes spanned by the tangents to such 
geodesies. The latter example gives insight to the geometry of the full manifolds of landmarks. 

1 Introduction 

In the past few years there has been a growing interest, in diverse scientific communities, in modeling 
shape spaces as Riemannian manifolds. The study of shapes and their similarities is in fact central 
in computer vision and related fields (e.g. for object recognition, target detection and tracking, 
classification of biometric data, and automated medical diagnostics), in that it allows one to recognize 
and classify objects from their representation. In particular, a distance function between shapes 
should express the meaning of similarity between them for the application that one has in mind, 
and also be mathematically sound and treatable. 

Among the several ways of endowing a shape manifold with a Riemannian structure (see, for 
example, [13, 14, 16, 21, 24, 25]), one that has recently gained popularity is inducing it through the 
action on the manifold itself of an infinite-dimensional Lie group of diffeomorphisms with a given 
metric, as described in [23, 26]. This approach can be used to provide a metric to several deformation- 
related shape spaces, such as the manifolds of curves [8, 22], surfaces [28], scalar images [2], vector 
fields [4], diffusion tensor images [3], measures [7, 9], and labeled landmarks (or "feature points") [10, 
11]. The actual geometry of these Riemannian manifolds has remained almost completely unknown 
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until very recently, when certain fundamental questions about their curvature have finally started 
being addressed [21, 22, 27]. 

This paper deals with the problem of computing sectional curvature for landmark points, and is 
based on results from the thesis of the first author [19]. The manifold of landmarks in Euclidean 
space is defined as 



This is among the simplest shape manifolds in that it is finite-dimensional, albeit with high dimension 
n = ND, where N is the number of landmarks, D is the dimension of the ambient space in which 
they live (e.g. D = 2 for the plane or the sphere). Therefore its metric tensor may be written, 
in any set of coordinates, as a finite-dimensional matrix. In fact it turns out that the inverse of 
the matrix defining the metric induced by the group action of diffeomorphisms (the cometric), has 
a relatively simple structure since each of its elements depends only on at most 2D of the ND 
coordinates. Hence the matrices obtained by taking first and second partial derivatives of the 
cometric have a very sparse structure — that is, most of their entries are zero. This suggests that 
for the purpose of calculating curvature (rather than following the "classical" path of computing 
first and second partial derivatives of the metric tensor itself, the Christoffel symbols, et cetera) it 
would be convenient to write sectional curvature in terms of the inverse of the metric tensor and 
its derivatives. So we have solved the highly non-trivial problem of developing a formula (that we 
call "Mario's formula" ) precisely for this purpose: for a given pair of cotangent vectors this formula 
c^xprcisses the corresponding sectional curvature as a function of the cometric and its first and second 
partial derivatives except for one term which requires the metric (but not its derivatives). 

The paper is organized as follows. We first give a few more details about the motivational 
example provided by the manifold of landmarks, and describe the metric induced by the action of 
the Lie group of diffeomorphisms. We then give a proof for the general formula expressing sectional 
curvature in terms of the cometric. This formula is used in the following section to compute the 
sectional curvature for the manifold of labeled landmarks. In the last section, we analyze the case of 
geodesies on which only two points have non-zero momenta and the sectional curvatures of 2-planes 
made up of the tangents to such geodesies. In this case, both the geodesies and the curvature are 
much simpler and give insight into the geometry of the full landmark space. 

2 Riemannian Manifolds of Landmarks 

In this section we briefly summarize how the shape space of landmarks can be given the structure of 
a Riemannian manifold. We refer the reader to [23, 26] for the general framework on how to endow 
generic shape manifolds with a Riemannian metric via the action of Lie groups of diffeomorphisms. 

2.1 Mathematical preliminciries 

We will first define a distance function : £^(K^) x £^(R^) R+ on landmark space which will 
then turn out to be the geodesic distance with respect to a Riemannian metric. Let Q be the set of 
differentiable landmark paths, that is: 



Following [26, Chapters 12 and 13], a Hilbert space (V, { , )v) of vector fields on Euclidean space 
(which we consider as functions — >■ K^) is said to be admissible if (i) V is continuously embedded 




Q:={q = {q\---,q 



,N 



) : [0,1] ^£^(K^) g« e Ci([0,l],M^),a 



i,...,n}. 
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in the space of C'^-mappings on M.^ which are bounded together with their derivatives, (ii) V 

is large enough: For any positive integer M, if cci, . . . , xm & and ai, . . . , au G are such 
that, for all ueV, {<^a, u{xa))^D = 0, then ai = . . . = um = 0- 

Thus (V, ( , )v) admits a reproducing kernel: For each a,x £ there exists £ V with 
(a, /(x))rc = {K^, f)v for all f eV. Further, (/3, K^{y))mo = {K^,K^)v which is a bilinear form 
in {a,l3) € (R-^)^, thus given hy a, D x D matrix K{x,y). 

In this paper we shall assume that K{x, y) is a multiple of the identity and is translation invariant: 
we then write K{x,y) simply as K{x — y)lD (where Id is the D x D identity matrix). There are 
other very natural admissible norms on vector fields v whose kernels are not multiples of the identity, 
e.g. one can add a multiple of div(u)^ to any norm and then K will intertwine different components 
of V. The most natural examples of the norms we will consider are given by inner products 

{u,v)v = {u,v)l := {Lu{x),v{x))^ dx, (1) 

where L is a self-adjoint elliptic scalar differential operator of order greater than D + 2 with constant 
coefficients which is applied separately to each of the scalar components of the vector field u = 
{u^ , . . . ,u^). By the Sobolev embedding theorem then V consists of C^-functions on which 
are bounded together with their derivatives. If ii' is a scalar fundamental solution (or Green's 
function [6]) so that L{K){x) — S{x), then the reproducing kernel is given by = K{ —x)a. A 
possible choice of the operator is L = (1 — A^A)'^ (where A <E M is a scaling factor, fc S N and A is 
the Laplacian operator), with fc > y + 1, in which case (1) becomes the Sobolev norm: 

•^^ e=l m=0 ^ ^ \a\=m 

When L = (1 — A'^A)'' the scalar kernel K has the form K{x — y) = 7(||a; — ?/||md), with: 

^(^)^ 2^.f-i.fF(fc)A^ (i)''^^--^^i^' 

where (with i/ = fc — ^) is a modified Bessel function [1] of order v. 

Now fix any admissible Hilbert space of vector fields. The space ^^([0, 1], V^) is the set of func- 
tions V : [0, 1] — >■ y such that: 



The space L'^{[0, 1], V) is a subset of L^{[0, 1], V) and is in fact a Hilbert space with inner product 
{u, 'y)L2([o,i],y) ■— Jq (W; ''^)v dt. It is well known [5] from the theory of ordinary differential equations 
that for any v G i^^([0, 1], V^), the D-dimensional non-autonomous dynamical system z = Vt{z), with 
initial condition z(to) = has a unique solution of the type z{t) = iJj{t,tQ,x). Let (p^^^{x) := 
ip{t, s, x); fixing t = 1 and s = we get := (p^i, which is the diffeomorphism generated by v. For 
an admissible Hilbert space we will call the set 

Gv := {^^ :veL'{[0,l],V)} 

the group of diffeomorphisms generated by V; by [26, Chapter 12] it is a metric space and a topo- 
logical group. But, in the language of manifolds, Qy is not an infinite-dimensional Lie group [15]. 
V is not a Lie algebra, but is the completion of the Lie algebra of C°° -vector fields with compact 
support with respect to || 
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2.2 Definition of the distance function 



For velocity vector fields v G -^^([0, 1], V) and landmark trajectories q G Q define the energy 



(4) 



/•I " J a 2 

E[v,q]--=l (hit, )||v + ^El^W-^(*'9"W)|L)^*- 

We claim that a distance function d on £^(R^) between two landmark sets (or shapes) / = 
{x^,x'^, . . . , x^) and I' = {y^, y^,. . . , y^) can be defined as 

d{I,I') := inf : v G L\[0,1],V), q G Q with q{0) = I, q{l) = /'}; (5) 

in the next subsection we will argue that the above function is in fact a geodesic distance with respect 
to a Riemannian metric. We treat the minimization of (4) as our starting point; it is the "energy of 
a metamorphosis" as formulated in [26, Chapter 13]. 

The above infimum is computed over all difFerentiable landmark paths q G Q that satisfy the 
boundary conditions (g"(0) = and = y", i = 1,. . . ,N), and vector fields v G L'^{[0,1],V). 

The resulting landmark trajectories {q°'{t),t G [OA]}a=i,....N follow the minimizing velocity field 
more or less exactly, depending on the value of the smoothing parameter X G (0,oo]; it is a weight 
between the first term, that measures the smoothness of the vector field that generates the dif- 
feomorphism, and the second term, that measures how closely the landmark trajectories actually 
follow the vector field. When A = oo we have exact matching, i.e. the landmark trajectories exactly 
satisfy the ordinary differential equations = v{t, g"), a = 1, . . . , N that arc obtained by setting 
the integrands of the second term in the right-hand side of (4) equal to zero. In fact, the exact 
matching problem can be equivalently expressed as the minimization of 



Eo[v]:= / \\v{t, Wydt 
Jo 

among all v e L'^{[0,1],V) such that '^^{x"-) = y°- , a ^ 1, . . . , N . When A < oo in (4) we have 
regularized matching, i.e. the landmark trajectories "almost" satisfy such set of ordinary differential 
equations; this allows for the time varying vector field to be smoother. For this reason the second 
term in (4) is often referred to as smoothing term; by allowing smoother vector fields the distance d 
is made tolerant to small diffeomorphisms and therefore more robust to object variations due to 
noise in the data. 



2.3 Minimizing velocity fields and Riemannian formulation 

By manipulating expression (4) we will now show that it is equivalent to the energy of a path q G Q 
with respect to a Riemannian metric tensor. 

Notation. Consider a landmark q = (q^,...,q^) in C^(R^). The D scalar components in Eu- 
clidean coordinates of the TV landmark trajectories = {q""^ , . . . , q""^), a ~ 1, . . . , N can be ordered 
either into a,n N x D matrix or in a tall concatenated column vector. We shall always use indices 
a,b,c,... G {!,..., N} as landmark indices, and i,j,k,... G {!,..., D} as space coordinates in M^. 
We will associate to each of the N landmarks q°- G R-'^^^ a momentum Pa G 'M}^^ (defined in the 
next Proposition) which we will write, in coordinates, as Pa = {Pai, ■ ■ ■ iPao), for each a = 1, . . . ,N . 
The components of momenta can also be ordered into an N x D matrix or in a long row vector. We 
chose superscript indices for landmark coordinates and subscript indices for momenta. 
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For a given set of landmarks (q^, . . . , q^) € wc will define the symmetric N x N matrix 

K(g) := (K{q°- — q''))^ (,^-^ ^- The matrix K.{q) is positive definite, whence invertible. 



Proposition 1. For a fixed landmark path q = {g" : [0, 1] — !■ M^}^ € Q there exists a unique 
minimizer with respect to v £ -^^^([0, 1], 1^) of the energy E[v,q], namely: 



N 



'{t,x) :=J2pa{t)K{x-r{t)), <e [0,1], xeR'', 



(6) 



a=l 

where the components of the momenta are given by: 

N 



Paiit) = ^ {K{q{t)) + ^f)~'^ ■ j/\t), t e [0, 1], 



(7) 



5=1 



a= 1,. . . ,N, i = !,...,£) (here Ijv indicates the N x N identity matrix). 



Remark. What the above proposition essentially says is that the vector field of minimum energy 
that transports the N landmarks along fixed trajectories is, at any point of time, the linear combi- 
nation of N lumps of velocity, each centered at a landmark point. The directions and amplitudes of 
the summands are determined precisely by the momenta. 

Proof of Proposition 1. Using property (ii) of the admissible Hilbert space V , [26, Lemma 9.5] shows 
that for given q= [q^,. . . ,q^) S £^(]R^) we have the orthogonal decomposition 

V = {v&V:v{q'') = 0,a=l,...N}®{v = Y.a=i(^aK{ : a„ e K^}. (8) 

Thus the minimizer must have the form 



N 



v{t,x) = J2(^a{t)K{x-r{t)), [0,1], X& 



(9) 



for some coefficients aa G C{[0, 1],M^), a = 1, . . . , A'", to be computed. For velocities of the type (9) 
the energy (4) can be rewritten as 



D N 



E[v,q] = / E E {^aiK{r - (i')au + A \a,iK{r - q') - q"\^] dt. 

^ i=l a,b=l 



(10) 



Setting the first variation of (10) with respect to coefficients Uai to zero yields the momenta (7). □ 
It is convenient, at this point, to introduce the ND x ND, block-diagonal matrix 



9{q) :- 







{K{q) + ^f)-' 



\ 



... (K(,) + ^)-V 



(11) 
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where the N x N block {'K{q) + ^y-) ^ is repeated D times; the choice of symbol g is justified by 
the fact that (11) is, as we shall see soon, precisely the Riemannian metric tensor with which we are 
endowing the manifold of landmarks, written in coordinates. 

Thus for a fixed path q € Q the minimizer of E[v,cj\ with respect to w € i^([0, 1],T^) is given 
by (6); since it depends on q we will write it, with an abuse of notation, as v*{q). We can define 

E[q]:=E[v*{q),q], (12) 

which depends only on the arbitrary path q G Q. The energy (12) is "equivalent" to the en- 
ergy E[v,q], in that: 

(a) if {v,q) minimizes E[v,q] then q minimizes E[q], and E[v,q\ — E[q\; 

(b) if q minimizes E[q] then {v*{q),q) minimizes E[v,q], and E[v*{q),q] = E[q]. 
Proposition 2. For an arbitrary landmark trajectory q G Q energy E[q] is given by: 

m= m^9m)mdt= ^^r(i)g'^(i)(K(g(t)) + f)^^di (13) 

^ ^ a, 6=1 i=l 

Proof. Following definition (12), formulae (7) for the momenta are inserted into the modified expres- 
sion (10) for energy E[v, q\. Simple matrix manipulations finally yield the right-hand side of (13). □ 

Remarks. Expression (13) has exactly the form of the energy of a path q with respect to Riemannian 
metric tensor (11). Whence given two landmark configurations / and /' in (MP) we have that 
if q minimizes (13) among all paths in (7 G Q such that g(0) = / and q{l) = I' then {E[q\)^/'^ is the 
geodesic distance between / and /'. By point (b) above we also have that {v*{q),q) is a minimum 
of energy E[v,q], so d{I,I') defined in (5) coincides with (i5[g])^/^ and is the geodesic distance 
between / and /' with respect to metric tensor g. 

The Lagrangian function that corresponds to energy (13) is: 

1 , N D H -1 

C{ciA) = \e9{cDi = \Y: T.ri''{^{<D + -f)^^- (14) 

a, 6—1 i—1 

In Hamiltonian mechanics [12] the "momenta" are defined as Pai = dC/dq""^, or, in vector notation, 
= dC/dq^^^ (for i = 1, . . . ,D). Applying such definition to (14) yields precisely equations (7) of 
Proposition 1. Whence the use of the term momenta is justified. 

Note that for small values of parameter A the metric tensor g, written in coordinates, gets close 
(up to a multiplicative constant) to the ND x ND identity matrix; in other words, for A — >■ 0, 
g converges to a Euclidean metric and the geodesic curves become straight lines. On the other 
hand, for A — )■ 00 (exact matching) the metric converges to [diag{K((7), . . . , K{q)}]^^ (block K{q) is 
repeated D times). In general, the block-diagonal form of the metric tensor g given by (11) follows 
from operator L in (2) being separately applied to each of the components of the velocity field; 
however the dynamics of the D dimensions of q are not decoupled since all ND components of q 
appear in each diagonal block of g. In the case of exact matching one can prove that landmarks 
"never collide" : in other words, it takes an infinite energy to make any two landmarks coincide. 

Figure 1 shows the qualitative behavior of geodesies in £^(M^), with A = 00. In the case 
illustrated on the left-hand side both landmarks travel in the same direction (from left to right, as 



6 



Figure 1: Two trajectories in £^(IR^). Bullets (•) and circles (o) are the initial and final sets of 
landmarks, respectively. The grids represents the two corresponding difFeomorphisms ipQi. 

indicated by the arrows): the two arcs of the geodesic "attract" each other, or in other words the 
two landmarks tend to "carpool" by using a velocity field with the smallest possible support so to 
minimize the part (i.e. the first term) of the Sobolev norm (2) of the velocity field. On the other 
hand when the two landmarks travel in opposite directions (as illustrated on the right-hand side of 
Figure 1) they try to avoid each other so that the higher order terms of the Sobolev norm are kept 
small; we shall return on the issue of obstacle avoidance at the end of this paper. A typical geodesic 
in £'*(E^) (again with A = oo) is shown in Figure 2. 

Conclusion. We have shown that distance d{I,I'), 1,1' E C^{M.^) defined in (5) is in fact the 
geodesic distance with respect to a Riemannian metric. In coordinates, the corresponding Rieman- 
nian metric tensor is given by (11), which is such that each element of its inverse (the cometric) 
depends on at most 2D of the ND coordinates. Whence the first and second partial derivatives of 
the cometric have a very sparse structure. This gives us motivation for deriving a general formula 
for computing sectional curvature in terms of the cometric and its derivatives in lieu of the metric 
and its derivatives, which will be done in the next section. 

3 Sectional Curvature in terms of the Cometric 

3.1 Generalities and notation on sectional curvature 

Let Ai be an n-dimensional Riemannian manifold. If we consider a local chart {U, ip) on the manifold 
with coordinates (x^, . . . , x"), we have the induced 1-forms dx^, . . . , dx" and coordinate vector fields 
{di :— gp:, . . . ,dn = afrr}- The metric tensor g : TM TM — ?> M can be represented as g\u — 
g{di, dj)dx'^®dx^ =: gijdx^®dx^ . We get a positive definite matrix with elements gij{p) = gp{di, dj). 
With an abuse of notation we will write gij(x) instead of [gij o ip^^)[x), x € fiU). 
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Figure 2: A typical geodesic trajectory in Bullets (•) and circles (o) are the initial and final 

sets of landmarks, respectively. The grid represents the corresponding diffeomorphism (p^i. 

Notation. We shall denote the partial derivatives of the elements of the metric tensor g as gij^k{x) := 
^gij{x) = dugij and gij.ktix) := gJl^k 9ij{x) = dud^g^y Also, we wiU indicate the cometric 
g~^\u — g^''di ® dj (so that g^-^gjk — S^) and their partial derivatives with g''-' f.{x) := ■T^g^^{x) 
and g"\t,i{x) := q^^9'^{x). 

For a tangent vectors X — X^di we consider the 1-form X^ := X^gijdx^ =: Xjdx^ (indices 
lowered), and for a 1-form a = aidx^ we have the tangent vector a" aig'^^dj (indices lifted). 

Indicating with X{M.) the space of smooth vector fields on the manifold A^, let V : X{M.) x 
A'(A^) — > X{M) be the Levi-Civita connection [12, 17] of the Riemannian manifold. The Christof- 
fel symbols F^^- are defined by Vg-dj — F*j(?fe, and it is well known that they have the form: 
rjj- = \g^^{git.j + gji.t — gij,i)- The Riemannian curvature endomorphism is the map R : X{M) x 
X{M) X X{M) -> X{M) given by R{X, Y)Z = \I x^yZ -^y^ xZ -SI \x,y\Z. In local coordinates 
R{di,dj)dk = Rfji^dt, and Rijkm '■= {Ridi,dj)dk,dm) = gmiRljk- The Riemannian curvature tensor 
acts on vector fields as follows: 

R{X, r, Z, W) = {R{X, Y)Z, W) (15) 

and in coordinates it is written as i? = Rijkmdx'^ ^ dx^ ® dx^ ® dx™'. The Riemannian curvature 
tensor has a number of symmetries: (i) RijM = -Rjiki', (ii) Rijki = -Rijtk', (iii) Rijki = RkUj] 
and (iv) Rijki + Rjkif. + Rkiji = (first Bianchi identity). With such conventions, the sectional 
curvature associated to a pair of non-parallel tangent vectors X and Y is computed by: 

= R{X,Y,Y,X) ^ R^.kmX^Y^Y'^X"^ 

'^^ ' > \\xp\\Y\\^-{x,Yy ||x|i2||r||2_(x,y)2- ^ °^ 

In order to express the numerator of sectional curvature (16) in terms of the elements of the 
cometric and its derivatives (i.e. g''^ , g^'' k ' ^^'^ S*"* ke) consider the covariant expression of the 
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Riemannian curvature tensor: 

which we call the dual Riemannian curvature tensor. Similarly we consider the covariant or dual 
Christoffel symbols 

■■= g'^^g'^QkuT^j , (18) 

which are symmetric in the indices r and s. 

To achieve notational compactness we will use the following symbols: 

Using that g = implies dkg = —Q~^ ■ dkQ ■ one immediately sees that 
Proposition 3. The following expression holds for the Riemannian curvature tensor: 

'^Rijkm — 9ik,jm H" 9jm,ik 9jk,im Qim^jk H" ik^ jm9oL^ jk^ imdot^ ' {^^) 

For a proof see [20, §24.9]. 
3.2 Mario's formula 

Proposition 4. The following expression holds for the dual Riemannian curvature tensor: 

^-K = - g - g + g + g + p >- a g p ^ a g 

I rX,u nv.s rX,u ns,v . u\,r iis,v uX,r iiv,s /'oi\ 

+ g g\p.g -g g\p.g +g g\p.g -g gxt^g (2i) 
+ g gxp. g + g gxn g - g gxn g - g gx/x g ■ 

Proof. We will manipulate (20) and write it in the form Rijkm = giugjrgksgmvR"^''^ by factor- 
ing giugjrgksgmv out of cach term; what will be left will be precisely the expression for R^^'^^ . 
The terms in (20) involving Christoffel symbols are, by (18): 

^l^^jM = gjrgmvg"^ g^^ = g^ug^rgk^g^v (r^r^ 5"") , (22) 

and similarly: ^fk^Lsc^i = g^ugJrgksgmv {T^^n-g"^) . (23) 

As we did in the proof of Lemma 18, ii g = Q^^ then djg = —Q~^ ■ djQ ■ and dmdjg = 
• {dmQ • • 9jQ + djQ ■ ■ dmQ — dmdjQ) ■ Q^^, i.e., in index notation, 

gik,jm = giu{g"^,m gxfi g^'^j + g^^j gx^ g^^^m ~ g^'',jm)gsk 

= g^ugksg,rg^. [5"«ff'"'(5"',, gx^ g^'\^ + 5"' « gxp ,9^% - g^^^^r,)] 

= giugjrgksgmv {g'^^'^gx,. g^''' + g^'^'^^gx^. g^''^ - s"''™) , (24) 
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where we have used definitions (19). Similarly, we can achieve the factorizations: 

gjm,ik = 9iu9jr9ks9mv [9 9\n 9+9 9\n 9 -9 j I (^5J 

-9jk,im — 9iu9jr9ks9mv [-9 9\n 9 " 5 9\ii 9 +9 )y V^Oj 

-9im,jk = 9iu9jr9ks9mv [-9 9XlJ. 9 9 9\n 9+9 )■ (^7j 

Inserting (22)^(27) into (20) we can write Rijkm = 9iu93T9ks9mvR^'^^'" , with i?"''™ given by (21). □ 
Proposition 5. The dual Riemannian curvature tensor may also be written as follows: 

2H =-9 -9 +9 +9 i^l) 

- ^{/^ g'^g''^, - 5".p (ff""'" + 5^"'") - (g"'''' + } (T2) 
+ 1 {5™,, g'^g^^a - {g'"''' + s"^'") - 5"% (s"'^'" + ^""'0 } (T3) 

- lig^""'' - s^^'OffA^ff""'" - 5'^"'") (T4) 

+ lig^""'" - s^^'OsA^ff""'' - 5'^^'") (T5) 

+ - g^^''')gUg''^'' - g'^''^)- (Tg) 



Proof. We will expand and recombine the terms in expression (21). The terms involving second 
derivatives need no manipulation and correspond to term Ti. The terms in the second line of (21) 
can be written as: 

5''^'"5A^5'"'''- 5"^'"5A^fl'*''''+5"^'''5A^fl''''''- 5"^'''5A^5''''''= (5^'"'"- 5^"''')5A^(5''"''- 5"''") 
which is precisely Tg. It is also the case that: 

p '-a 9 - 9 gxfig - g gxp.g'^ 
= 5 [(5^''" + 9^"'! - 5™'^ gxp g'^gai. [(5'^"'^ + g^'n - - ^''^'''sa^ 5''''" - 5"^''5a^ 5'^"''' 

+ Ua^'^''' + g'^'lgxM^'^ + fm ~ g'^^'^'gx^g'''''' - g'^^'^gx^.g'^'"''' 

= T3 + K^^"^'" - 5^'''")5A^(5'^"'' - 5"''") = T3 + T5. 

Similarly one can prove that: -2T'p^T^'' gP" + g'^^'^^A/i 5'*'''" + g^'^'^gx,^ g^'''' = T2 + T4. □ 

For an arbitrary pair of tangent vectors X — X^di, and Y = Y'^di in Tpjvi wc consider the 
covectors X^ = Xidx^ and Y^ = Yidx^ in T*M.^ with = gijX^ and = gijY^ . The numerator 
of sectional curvature (16) may be rewritten as RijkmX'-Y'^Y'' X"^ = R^'''''" XuYrYsX^. 

Theorem (Mario's formula). For an arbitrary pair of vectors X = X^di and Y = Y^di in TpM 
the numerator of sectional curvature (16) at point p G M. may be written as: 

g{R{X,Y)Y,X) = R^^'^ XuYrY.X, = 

= {X^Yr - YuXr) (iff^"'™ + 15"% fl'"^'" - i 5"V 5™'" - I 5^"''^5A^ g""'-") {XsY, - Y,X,). 
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Moreover, if we extend X and Y locally on Ai to constant 1-forms in terms of local coordinates 
(i.e. make its coefficients X„.Yr constant functions), then the formula becomes: 

g{R{X,Y)Y,X) = 

= [lXX{\\Y^f) + lYY{\\X^f)-l{XY + YX)g-\x\Y^)} 

+ {i||%-l(x^y^))f -i5-'(d(||x^f),d(||y^in)}-|ff([x,r],[x,y]), 

wliere the term, in the first set of braces equals the sum, of the first two t,erm,s in the coordinate form, 
the term in the second set of braces equals the third term in the coordinate form and finally the last 
terms are equal. 

Proof. We will write the six terms provided by Proposition 5 as T"*"*^, i = 1, . . . , 6. We have: 

Tursv v V V V „us,rv y V V V rv,us y V V y i ,,rs,uv y V V y i ,,uv,rs y V V y 

= 9"'''^^ {^X^YrYgXy — XpYuYyXs + XrYuYgXy + XuY^YyXg) 
= - YuXr){XsY, - YsXy), 

where the second step follows from relabeling the indices. As far as T2 and T3 are concerned, 

{TY'-+TYnXuYrY,Xy = gf^g^^.X^X, - F.X.g™^ gP-g^^^^X^Y,} 

+ l{YrY, g^^pigP^'^ + + X^X, g^^^ig"''' + ff^^'^l^n 

- YrX,. g^\p{g'^^' + gP'''')XuYs - X^Y, + g'"'''')YrXy] 

= - i{2i;n.9'"% gP'^g^^.XuX,, - 2y,X,g™^^ .g^-.g^^X^K,} 

+ \{2YrY, g^\^ g'^'^^XuX, + 2XuX, g-\, g'^'^'Y^Y, - IX^Ys + ^'"''O^rX^ 

= — I; 9^^\p gP" g^'' ,y\YyYsXrX^ + YrY^X^Xs — Y^XyX^Y^ — YyXrXgYu^ 

+ 6'"*,p g''^'" {YuYgXrXy + XyXgYrYy — XyYgYrXy — XgYyYyXr^ 

= ( - 3 3™'" + {XuYr - YuXr){XgYy - YgXy), 

where, once again, step (*) follows from relabeling the indices. Also, one can easily see that 
TT'^'XuYrYgYy = -lYrY{g^-^'^ - g^'^gx^''''' - g'^^'^^uXy = 0. Finally, 

(^T^ursv _^_r^ursv^X^YrYgYy 

= \ Y,Xy[g^-^- - g^^n9x,{cr^' - g^'^^uYs + Y,Xy,{g^^'^ - g^^ngxy.ig"'''' ' g^^'^XvYs 

= f y.x„(/'^'" - g^^^ngxM^'^' - g^nXvY 

= — I g^^'^ gxf^g^'^'" — YuX,i.XsYy + YyXj.XyYg + Yr-X^XgYy — YrX^XyYg^ 

= - i g'''''''gXf,g''''^XyYr - YyXr){XgYy - YgXy). 

Divide by 2 to get the coordinate formula. The non-local version of the formula follows easily by 
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bringing the X and Y's into the formula; thus: 

= XrX,{{\\Y'f),,,g>'^g''^ + {\\Y' f) g"- g'^^) 
= X,g'^^{Xrgn\y'f),p),„ = X^{Xf\\Y%)^^ = XX{\\Y'f). 
A typical term from the third part of Mario's formula is rewritten like this: 

the other terms arc similar. Finally, it is the case that: 

{X^Yr - YuXr)g^^^'-dx = (XuYr - YuXr)g^^^^ g^'-^A = {XuY'' - ^^n^"),?^" „ 9a 

= {{Xug^n^r, Y^ - {Y^g^"),r,X^)dx = {X\Y^ - Y\X^)dx = -[X,Y], 



and the proof is easily completed. □ 

Remeirk. It is convenient to split Mario's formula in four terms: 

i?i := i {XuYr - - Y,X,) , (29) 

R2 ■■= \ {XuYr - yuXr)g^\p .g^'-'" - Y,x:) , (30) 

i?3 := liXuYr - Y^Xr) ( - i 5"^. 5™'") {X,Y, - Y,X,), (31) 

i?4 := liXuYr - Y^Xr) ( - 1 ^^"■'•flA,. i?''^'") {X,Y, - Y,X,) ; (32) 



all the terms with the exception of R4 (where g appears, but not its derivatives) depend only on 
elements of the cometric and their derivatives. 

Remark. The denominator of sectional curvature (16) can also be expressed in terms of the comet- 
ric: 

||yf - = x„x,r,n(5"«5™ -5™5-). (33) 

4 Curvature of the Manifolds of Landmarks 

In this section we will apply Mario's formula to the computation of sectional curvature for the 
Riemannian manifold of landmarks, introduced in section 2. 

4.1 Hamiltonian formalism 

On the A''D-dimensional manifold C = £^(]R^) of landmarks wc consider the Riemannian metric g 
given, in coordinates, by the matrix (11); it is in block-diagonal form and we write its generic 
element as g(ai)(bj)^ with a, & = 1, . . . , (landmark labels) and i,j = 1, . . . , D (coordinate labels, 
respectively of landmarks a and b). More precisely: the matrix 17(0') is made of D square {N x A'') 
blocks; indices i,j = !,...,£> indicate the block, whereas indices a,h= 1,. . . ,N locate the element 
within the («, j)-block. Therefore if we indicate with hab{<l) the generic element of the N x N matrix 

(K(g) + ^)"^ we have that 

9iai){bj) = hab{q)Sij, a,b=l,...,N, i,j = 1,...,D, 
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where 5ij is Kronecker's delta. Similarly, if wc indicate as the elements of the cometric 

tensor g{q)-^ , they are given by .g(")(''^)((j) = (5*^ where = K{q''-q^) + ^. In analogy 

with the notation introduced in section 3 we also denote the partial derivatives by g^°'^^'^^^^ ^^^^ = 
and = g^Hf^S^"*^^"^'^; they will be computed later. 

For simplicity from now on we shall assume that A = oo (i.e. that we are dealing with exact 
matching of landmarks); so the element of the cometric becomes = K{(f' — q^)5'^^ . The 

Hamiltonian function [12] for the system can be written as: 

^ N D ^ N D 

a, 6=1 i,i=l a, 6=1 i,j=l 

that is n{p, g) = 5 E^6=i Kiq" - q^){pa,Pb)^D- 

Proposition 6. Hamilton's equations for the Riemannian manifold of landmarks are: 

N 



6=1 



^ a=l,...,N. (34) 

Pa = -Y.VK{q'^-q'){p,,Pb)^, 

6=1 



Proof Equation (7) can be written as = X^fe^i K{q'^ — q'') pu^ for a = 1, . . . ,N, i = 1,. . . ,D; 
alternatively, computing q'^^ = yields the same result. Also: 

gl^K{q^' - q^\ . . . , g*"^ - g'^^) = Ef=i ^{<i' - - <f') 

= Ef=i l^il' - in {S'a - m = ^{q" - g^) {St - 51) (35) 

so that 

Pai = -^{p,q) = -i Ef=l 1^(9" - Cl^) {Pa,Pc)^o + \ Ef=i if (g*" - g") (p6,Pa)M« 
= -Er=l|f(9"-9*)(Pa,f6)M-, 

where the last step follows from the skew-symmetry of V-ft'"'' in indices a, b. □ 



Corollary 7. IfPa{to) = for some landmark a = 1, . . . ,N and time to € K, then Pa{t) = 0. 



4.2 Notation 

Let K{x), X e be the scalar kernel that defines the metric; we assume that it is twice continuously 
differentiable and symmetric, K{x) = K{—x); for now we shall not assume rotational invariance. 
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We define: 



d,K{x) := —^{x), d.K-" := d,K{q- - q") G M, 

VK ~ {diK,--- ,dDKf, WK"'' :=WK{q'' -q^)&mP, (35) 

ZJ^iv- := Hessian {DIK), D'^K"^ := - q^) G IR^^^. 

Note that ^K'^'' = -V/f*", V-fC"" = and D'^K''^ = D'^K^'', for all o, 6 = 1, . . . , TV. 

For a fixed set of landmark points q in C = (UP) consider any pair of cotangent vectors 
a, /3 G T*L: we shall write a = (ai, . . . , qjat) and f3 = (/3i, . . . , /Jat), where each component is D- 
dimensional. We define the vector field a^°^ : IR^ — >■ and its values at the landmark points by: 

N N D N D 

a^^'{x) := K{x - q^)ab = ^ ^ ii'(a; - q^)abjdj = ^ ^ ^^(a; - q^)au5'^dj, x G K^, 

6=1 6=1 j=l 6=1 i,i=l 

AT 

(a")" := a^°''(g«) = ^ i^'^^'afe, 
6=1 

which are, by virtue of formula (6), the velocity field 0^°^ on induced by the landmark momen- 
tum a = (ai, . . . jUn) and the corresponding landmark velocity a" G Tq£ (which obviously coincides 
with the second of Hamilton's equations (34)). Note that a" = {Oi,. . . ,a^j^) is the tangent vector 
in Tg£ with metrically lifted indices. 

The curvature of the Riemannian manifold of landmarks will be expressed in terms of three 
auxiliary quantities which wc now introduce. Wc will call those force, discrete strain and landm,ark 
derivative. We start with the force. For a fixed covector a = (ai, . . . , ajv) € T*C, having the dual 
vector extended to a vector field a^°^ on all of MP allows us to take its derivatives at the landmark 
points, a D X D matrix- valued function on M.^: 

N 



{Da^^^y^ix) := diia'^^yix) = ^ ab,-5ii^(a: - g"), 

6=1 

N 

{Da'^°y{q^)=J2diK''''^bj. 



6=1 



For a trajectory {q{t),p{t)) of the cotangent fiow one has that {pi{t), . . . ,PN{t)) G T*^^^jC for all t 
whore the trajectory is defined, so the above notation can be used to rewrite Hamilton's equations 
in a more compact form. In particular, the following result holds. 

Proposition 8. The second of Hamilton's equations (34) can be written as 

Pa = -Dp''°\q'')-Pa a = l,...,iV. (37) 

Proof Pa,. = -Eb=ldiK''HPb,Pa)RO = - ^f^l {Eb=l ^iK'''' Pbj)Paj = ' Ilf^liDp'^niPaj = 

-{Dp^°'{q'')-Pa)^, for any a = 1,..., TV and i = !,...,£). □ 
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For a fixed cotangent vector a G T*C, this motivates defining the right-hand side of (37) to be force: 

Fa{a, a) := Da^'''{q'') ■ aa, a=l,...,N. 

The full bilinear, symmetrized force may be thought of as a map F : T* C x T* C ^ T* C. We call 
the covectors given by this the mixed force, with the definition: 

Fa{a, /3) := i (£»a^-((7«) • /3„ + £»/3^-(g«) • a„), 

Fai{a,l3) := '^J^T.^i^'^' ("""jf^aj + Phja^j) = i ^5iif'^''((a„,^b) + (/3a,a(,)), (38) 

j=l 6=1 6=1 

for a = 1, . . . , TV and i = 1, . . . ,D. (The angle brackets are inner products in M.^.) Note that the 
"complete" cotangent vectors a = {ai, . . . , a^r) and (3 = {f3i, . . . , (3^) (not only their a-components) 
are needed to compute each component Fa{a,(i) of the mixed force. The mixed force has simple 
interpretation. If we extend a and /3 to constant 1-forms on C, then the differential of the map 
q ^ g~'^ia,P) = Ea,6^(9" - <l^){<^a,Pb) is given by: 

AT D 

dg-\a,(3)= J2 J2^iK{q''-q'){dq'^'-dq'"){aa,l3b) 

a, 6—1 i—1 

N D 

= T.J2 Wq'' - q") ((^a, (3b) + {0a, dq''' = 2F{a, (3). (39) 

a, 6—1 i—1 

For a fixed a € T*C the second component of the curvature formula, called the discrete vector 
strain, is defined by: 

JV 

Sabia) := {a^r ' or S^biaY := J^i^'"' ' K'"')a,i 

for all a, 6 = \,...,N. These are vectors and are skew-symmetric in the points a,h: Sab{ce) = 
—Sha{oi), Saa{oi) = 0. The scalar quantities: 

JV D 

c=l i=l 

we define to be the scalar compressions felt by kernel K; they are symmetric (since both factors in 
the inner product are skew-symmetric), i.e. Cab (a) = Cba{oi), with the property Caaia) = 0. We 
call these compressions because if X is a monotone decreasing function of the distance from the 
origin (the most common case), then VKab points from q°- to q^. 

Finally, if v and w are any two vector fields on landmark space, we may write their Lie derivative 
as the difference of covariant derivatives: 

b,t«]z: = Vf'*i^*(«;)-V£'«^*(^;) 

where the flat connection on £ is just the one induced by its embedding in M.^^. In other words, 
y£,flat^^-j jg ^jjg ^g^al derivative of w in the direction v if we use the coordinates on landmark 
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space: that is, V^'«^*(w) := Eai ^("'"O^ai = J2aiJ2bj i9bjw''')dai. If a, p are constant 1-forms 
everywhere on we can take v = and w = now as vector fields on C, and then we find: 

= E E/"")" ^J-^" (^fc" - ^ft') dai = Y.Y. - Pbi dai 

a,i b,c,,j a^i b^j 

= E E - '^K''}^^b^ ^a. = E ( E ^a6(a)/3fcz)a.i. 

This is a vector in TqC which we define to be the landmark derivative of /J" with respect to aK The 
coefficients with respect to dai , dao (for fixed a) are the elements of the following vector: 

AT N 

^''(a,/?) :=ECa6(a)^fa= E(^'''-^'")<"<='^^"')^'" a=l,...,N. (40) 

6=1 6,c=l 

We have that D{a, P) = {D^ia, is the A^D-dimensional vector of the coefficients of Vfi^''*(;3«) 

with respect to the basis {dai} of T^C. In particular, the coefficients of the Lie bracket of a" and /3' 
as vector fields on C are given by D(a, (3) — D{/3, a). 

4.3 General formula for the sectional curvature of £^(M^) 

We can write sectional curvature of £^ (R^) in the following way, where we have split it in the terms 
introduced by (29)-(32). From now on ( , ) will indicate the dot product in M^, while { , ) 
and ( , )t'C will be the inner products in the tangent and cotangent bundles of £ = £^(M^), 
respectively. 

Theorem 9. The numerator of sectional curvature o/£^(R^), for an arbitrary pair of cotangent 
vectors a and (3, is given by R{a^ , (3^ , , a^) = E^=i 'with: 

^1 = 5 E ® -13a® Sab{a)f{lD ® D^K'^') {ab <E> Sab{P) -Pb® Sabia)), (41) 

a^^b 

i?2 = E ((^"(". ^a(^> P)) + (D^W, P),Fa{a, a)) - (D«(a, (3) + D-{p, a), Fa{a, p))) , (42) 

a 

= E {{Fa{a, 13), F,{a, /3)) - (F,(a, a), F,{/3, /3))) , (43) 

ac 

Ra = -l\\[aKpHc\\lc = -|||£'(a,/3) - £'(/3,a)||J,_,. (44) 

In the formula for the first term i?i we have used the well known definition of tensor product: 
(vi (g) V2)'^{Mi (8) M2){wi (g) ■W2) {vf MiWi)(vJ M2W2) , and in the formula for the fourth term R4 

the norm for D x matrices || J||^ := Y.f=i J2a,b=i J^aJ2bAab- 

The theorem is proven by applying Mario's formula to the cometric of the manifolds of landmarks. 
One needs to compute the elements of the cometric and its derivatives in terms of the kernel and 
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its derivatives (36). In agreement with notation (19) we will define (note that we will keep using 
Einstein's summation convention wherever possible): 

{ai){bj),{de) (ai)(bj) {ck){M) J (ai){b3),{ck){de) (ai)(6j) {ij.p){ck) {ia){M) 

^ ' ^ tick) y y • y ,ii^p){$<^) ^ ^ 

Lemma 10. It is the case that 

5^"'*'"l(cfe) =dkK-'' iS^,-St)S^\ (45) 
5^"^^^'"l(cfc)(..) = dleK'^' i^c - S'J (<5S - 5^) <5^^ (46) 

g{ai){bj),ick){dl) _ Q'^^J^^b ^J^ac _ J^bc^ (_K'"<^ — K*"^) 5'-'. (48) 

Proof. Since gi"-^)^) = K°-^5'^ and also ^K{q'' - q'>) = dkK<'''{S^ - S^) by (35), equation (45) 
follows immediately. Similarly to (35) one can prove that g~adkK{q'^ — q^) = dj^.K"''' {6°i — 6^), 

whence: = gf^ 9^"'^^''\^,k) = ^IkK"' (^3 - ^^)Vc" - ^c') so (46) holds too. Now, 

by expression (45): 

g(ai){bjUdl) ^ g{ck){de) ^ Q^j^ab (^a _ ^6) ^ij j^cd §ke ^ Q^j^ab (^j^ad _ j^bd^ §ij _ 

which is (47). We can use (46) to compute ^(ai)(fej),(cfe)(d«) = g^"''^'^^^\^^p)(^^^^ g{i^p)(ck)g{ia)(de). 
g{ai)ibj),ick)ide) ^ Y.f,p^a dpaK"'' {6^ - (5|^) ((5| - <5|) S'^ K'"' Sp'' K^'^ 5"^ 

which completes the proof. □ 

Proof of Theorem 9. We will compute terms R\,...,Ri introduced by formulae (29)-=-(32). For 

simpUcity, sometimes we will write Da^^ instead of D(x^°'^{(f). 
• Computation of Ri. We have Ri = ^{a « - ^ -(""^^^^^ 

expression (48) into such formula yields: 

2^1 = Eall indices («a«/3cr - Paua^r) d^^K^'' {K'^' - K^^) {K'^'^ - K^") (5"« {a^s^dv ~ Pbs^dv) ■ 

Performing the above multiplications gives rise to four terms, which we will now compute one by 
one. First of all we have: 

2^1,1 := Eall indices c^auficrahsPdv d^.R-^ (if - K^-^) (K^^ - K^^) S"' 

= Eafc^r^fcE™ [Sabm^dXK-b[Sam]^ = Eab^I^b {sMfD^K-^sM\ 

= Eab ("a ® Sab{P)f{lD ® K^") {a^ ® Sab{P)) I 



• Computation of Ri. We have i?i = ^{aaul3cr - PauCtcr) 9^°''"^^ ^^'^''''^^'^'"\absl3dv - l^bsCtdv)- Inserting 
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similarly, 

2i?i,2 := - Eaii aauPcrPbsadv d^^K'^'' (K'^^ - K'"^) {K^'^ - K^'^) 6^' 

2i?l,3 := - Eall PauacrabsPdv O^^K'^" (if - K"-) {K'^" - K^") 5^' 

= -Eab (A. ® Sab{a)f{lD <E> D^K-''){ab SabiP)), 
2i?i,4 := Eall PauC-^crPbsadv d^K-'' (K-- - K^-) [K-^ - K"'') 

= Ea5 (/5a ® Sab{a))'^{lD ® K'^'') (/3fc ® Sab{a)). 

Now we can take the summation Ri = Ei=i which yields precisely expression (41). 
• Computation of R2. We may combine equations (45) and (47) from Lemma 10 to get; 

= drK"^ [d^K'"'{K'"^ - K""^) - dyK^'={K^'^ - K'='')] S'^". (49) 
Inserting (49) into 2i?2 = (aauPcr - /3a«acr) 9^^''^^""^'^''"\c,bsPdv - Pbs^dv) yields: 

2^2 = Eall indices { ^auPcrabsfidv OrK^" [d^K^^K^" - K^") - d^K'^K"" - K^')] <5«^ 

- f3auaMf3dv drK-'[d,K-'{K-^ - if^^) - d^K^'^K''^ - if^^)]<5"^ 
+ (3auaer(3bsadv drK^' [d.K^^K'^'' - K-'') - d,K'>-{K'"' - K-'')\5^'], 
which immediately implies: 



i?2 = 



\ Eabcdi^a, ab){Pc, ^K-^) [{Pd, VK-%K»<^ 




{0d,VK^''){K'"^- 




{.=■■ 


^2,1) 


-h^abcdi^a, 13b) {13c, VK-^) [{aa, VK-'^){K-'^ 




{ad,WK^%K^'^- 




(=: 


-^2,2) 


- \ Eabcd (Pa , a,) (a,, Vif [(/3d, Vif-) {K^^ 




(/3d,VX^^)(if^'* 




(=: 


^2,3) 


H Eabcdil^a, ^b){ao, VK-') [{ad, WK-''){K-^ 




(ad,VJs:'"=)(Js:'"^ 




(=: 


-^2,4) 



We will now manipulate terms R'2,i,- ■ ■ ,R2,4 one by one. Since V-fC" = — V.ft' , by relabeling the 
indices wc have 

^2,1 = Y.abcdi^a,ab){l3e,VK-'>){Pd,VK--){K-'' - K'-i) 

= J:abMa,ab){Po,VK-»){j:,K-'^p, - Y.dK'^Pd,^K--) 

= T.abMa,ab){Pc,^K-''){SM,^K-'^) = Eabc{»a,ab){Pc,'^K-b)CM 
= Eab{^a,ab){EcCac{P)Pc,'7K-^) = EabK, ",)(Z)«(A /3), Vif"^) 

= Ea6 D"{P, P^VK-'alaa = Ea D-{p, pYDal-^ ■ a„ = Ea(^"(/3. P), Fa{a, a)). 
Similarly, i?2,4 = '^^{D"{a,a), Fa{(3, (3)) ■ It is also the case that 

^^2,2 i Ea5c("a, /3b){l3c, Vif [(Ed (^"' " V K^-) - {J:^{K'"' - K-<')ad, VK^")] 
= -5 EabMa, Pb){l3c, VK-") [{Sac{a), V - {Sb,{a), VK^'')] 
= - \ EabMa, Pb}{f3c, VK-"} [Cacia) - C^] ; 
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relabeling the indices (and using the fact that VK""^ = —VK^"") yields: 

^2,2 = - \ Eatc [i^a, Pb) + Pa)] {Pc, VK-^)C,,{a) 

= - 5 Eafc [("a, Pb) + {at, Pa)] (EcCac(a)/3e, Vi^«*) 
= - 5 Eafa [("a, Pb) + {ah, Pa)] {D'^i^, P) ,V K^^) 

= - \ ZabD^ia, P)'^ [VK-^Plaa + VK-'alpa] 

= - \ Ea^"(«> {^Pl°" ■ "a + J^a'r ■ Pa] = -Y.a{D''{a, P),F,{a, /?)). 

Similarly, i?2,3 = -J2aiD''{P,a),Fa{P,a)). By the symmetry of 

R2,2 + i?2,3 = -Ea(^"("> P) + D^'iP, a), Fa{a, fi)) 

Adding the above sum to the expressions for ii2,i and i?2,4 finally yields (42). 
• Computation of R3 . We have 



But by Lemma 10 

,(?7cr) y Z^r] 



= {'VK"'', Vii''='^)^"^5™(if"'= - K'"^ - K^" - K*"^), 

whence: 

-8i?3 ^ Eaii { (Vii^''^ VK'"^){K'"' - K'"' - K^" - K'"^) 

■ {oiauPcrabsPdvS'^^S^^ -aauPcrOtbsPdyS^'^S'''' -aauPcrOtbsPdv^^^'S'''' + aauPcrOibsPdvS"''S'~'')} 
= Ea6cd{(("a'"b)(^c,/3d) - {O-a, Pb){Pc,OLd) - {Pa, ab){ac, Pd) + {Pa, Pb){ac, (Xd)) 

Relabeling the indices in the above expression yields: 

-8i?3 = J2abcd [^{aa,ab){Pc,Pd) " ^{aa, Pb) {Pc, ^d) " ^Pa, Otb) {oLc, Pd) 

- 2{aa, Pb){Pd, ac) - 2{l3a, ab){ad, Po)] {VK''', VK^'')K-'' 
= EabcdK'"'[Salab{VK->>)^WK'^''PJPc - 2alPb{VK-^)'^S/K^''alP, - 2plab{VK-^rVK^''Pla, 

- 2alPb{VK'^YVK^''pJa, - 2l3^ab{VK'^YVK'=''aJp,] 

= EacK'"'MiDa^°YiDp'^nPc - 2aliD(3^'>YiDa^nPc - 2(3^{Da^°Y{Dp^"^)ac 

- 2«^(i^;9r)^(i?/3^)a, - 2/3j(i^ah-)^(i?a^-)/3,] 

= Eac [^{Da'r ■ aa, Dfi'r -Pc,)- 2{Da'r ■ Pa + Dpl°^ ■ aa, Da"^-^ ■ + Dp^°'^ ■ a,)] . 
= Eac K''' [&{Fa{a, a),Fc{P, P)) - 8(Fa(a, P),F,{a, /?))] , 

which is precisely (44). Alternatively, this can be derived from formula (39). 
• Computation of R4. It is the case that: 
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By Lemma 10: 

- /3a« [E.5rif«"(Ec^*'«cr)] + Pan [Er^ri^^" (Ec^"'"-)] } -^^^ 

= Ea„ {«a« [(Vi^«", - (Vif ^^ - /3a« [(V^^«^ a^-) - (VX«^ } 

So if we define tlie matrix Hia := J2b ['^ab{^)ctbi — Cab{c()Pbi] , i = 1, ■ ■ ■ , D, a = 1, . . . , N we have: 

o „ n AT „ 

«s {r; Aju u=l ^,?7=1 

Alternatively, this can be derived from formula (40). □ 

The denominator (33) of sectional curvature for £^(]R^) is given by the simple formula: 
Proposition 11. For any pair of cotangent vectors a,l3 G T*C, 

M\%'C WWl'C - («./3)t-£ = E^''^''(("-"f)(/3c,/3d) - (aa,/36)(ac,/3d)). (51) 

abed 

Proof. Using double-index notation we may write equation (33) as follows: 

= Eabcd /3cr Pdv {K^H^^ K^H^^ - R'^d^uv j^bcgsr^ 

and (51) follows by relabeling the indices. □ 

4.4 The rotationally invariant case 

Finally, suppose the Green's function K is rotationally invariant: 

K{x) =-f{\\x\\), X gR^, with7e C2([0,oo)). (52) 

Whenever this holds, we will use the convenient notation: 70 := 7(0), 'jab '■= 7(||5° — <l''\\), I'ab '■— 
- ^ll), and 7^'b := i'{\\q^ - q^\\) for a, 6 = 1, . . . , AT. Then we can evaluate the first and 
second derivatives of K: 

Lemma 12. For rotationally invariant kernels, it is the case that: 

VK{x)=i{\\x\\)^^, (53) 



x\ 



D^Kix) = 7"(l|x||) - ^1 ^ + « 1. (54) 



= 7"(INI);^i + ^Pr-(.), 
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where Id is the DxD identity matrix and Pr"'"(x) := Id — is projection to the hyperplane o/M-^ 
normal to x. 

Proof. We have that diK{x) = 7'(||2^||) and (53) follows immediately. Also, 
= ^gf,y(IN|)+7'(INI)^gfia.^+7'(INiy^ 

= 7"(l|x||)f^ + ^<5- - y(l|x|i)f^ = [7"(||x||) - ^] + «<5-, 
which implies (54). □ 

Because of (53), in the rotationally invariant case, the "scalar compression" Cab{oi) really does 
measure a multiple compression of the flow a" between and q^. We can decompose the vector 
strain Sab{oi) into the part parallel to the vector g" — and the part perpendicular to this: let 

Uab ■■= ||gll^6|| and define 



Slb{<^) ■■= {Sab{a),uab), S^,{a) := Sab{a) - Sl{a) u^b- (55) 

Note that S^^^^^a) is a 5ca/ar while S^f,{a) is a vector. In particular we have that Cab{a) = I'ab'^abi^)- 
Moreover, formula (54) allows us to simplify the first term Ri in the curvature formula. Substituting 
(54) into (41), we get the rotationally invariant case for Ri: 

Proposition 13. In the rotationally invariant case (52), we have that 

^1 =E - ^i(/3)"a> sl{a)Pb - sl{P)ab) (56) 

^^^%^(5^,(a) ^(3a- S^M ® aa, S^.ia) Sj^,{P) ® ab)) . 



ai^b 



Proof. For any pair of covectors r] and in T*£, by (54) we have that: 

Sab{Tl) D^K'^" SM = 7;'b Sabilf Uabul, Sab{fl) + ^abilf Pr^(«a6) Sab{fi) 

= <b Slb^ri) + p^;^ {SiM, 5i,(M)). 

Inserting this expressions into (41) yields the desired result. □ 
4.5 One landmctrk with nonzero momenta 

A simple special case is when only one landmark carries momentum. We now compute the numerator 
of sectional curvature when both cotangent vectors are nonzero at only one of the £>- dimensional 
landmarks (g^, . . . , g^). We define: 

(T;/:)i := {q e TIL | = for a > l} 

so that the elements of the above set are cotangent vectors of the type r/ = (?7i, 0, . . . , 0). 
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-i) — (I — (I — (I — (I — (I — () — () — <>- 

-i lVu <) () () () <) () ( >- 



^1 — CI (I — (I — (I — (I — (I — (I — <K 
HI — (I <l — (I — (I — (I — (I — (I — ih^ 



HI (I <l — (I — (I — (I — (I (I — (1^ 
HI — CI — (I — (I — (I — (I — (I — <b,4K 



HI — CI (I — (I — (I — (1^ 



HI — CI — CI — CI — CI — CK 



HI — CI — CI — CI — CK 
HI — CI — CI — CI — CI — C» I I O I O 



HI — CI — CI — CI — CI — CI — CI — CI — CK 
HI — CI — CI — CI — CI — CI — CI — CI — CK 



HI — CI CI — CI — CI — CI — CI — CI — C>^ 



-<l — (I — CI — C>- 
-C l I H I C H-<>- 



^1 — C) CI — CK 
^1 — CI it ih^ 



— C> CI Cl^ 

^l — CI — CI — CK 



HI — C> CI — CK 



HI — CI — CK 

^l — c> CI — 



^l — CI — CI — CK 
^l — CI — CI — CK 



o o o o 



o o o ° 




Figure 3: Dragging effect of one momentum-carrying landmark (bullet •) on a grid of landmarks 
(circles o), with 7(x) = exp(— i^), a = 1.5. Left: initial configuration, with initial momentum 
Pi — (2.7,1.8) also shown. Right: configuration after one unit of time, with trajectory of also 
shown; the red grid represents the diffeomorphism ipQi, obtained by integrating a^™ in time. 



Proposition 14. In C^{R^), for any pair a, /S e {T*C)i the four terms o/ i?(a'*, /S", a") are 
given by Ri ^ R2 ^ R3 ^ Q and R^ = -| J2a,b=2(^a, Hb)jg,D{K-^)ab, where 

Ha (7ai -7o)((ai,Vif°i)/3i - (/3i, Vif'^i)ai) , for a > 1. 

Proof. Using formula (38), we see that all mixed forces are zero. Therefore, the result for 
the first three terms follows from Theorem 9. Also, by (40), D''{a,(i) = (7^1 - 7o)(77i, VK°i)/3i 
since a, /3 € {T*C)i\ a similar expression holds for D°-{f3, a), which concludes the proof. □ 

Therefore when a,f3 € (T*C)i the sectional curvature is always negative; we can understand 
this by considering the geodesic flow in this case. It follows immediately from Proposition 6 that if 
we start with zero momenta pa at all g", a > 1, then the momenta at these points stay zero, while 
the momentum at q^ remains constant. Thus the velocity of q^ is just given by K{0)pi and this is 
constant. The point q^ carrying the momentum moves in a straight line at constant speed, while 
the other points q"^ (a > 1) are carried along by the global flow that the motion of q^ causes and 
move at speeds = K°'^pi, which are parallel to q^ (but not constant). As shown in Figure 3 (the 
central landmark q^ is the only one carrying momentum) what happens is that all other landmark 
points are dragged along by q^, more strongly when close, less when far away. Points directly in 
front of the path of q^ pile up and points behind space out. 

Negative curvature can be seen by the divergence of geodesies. If you imagine slightly changing 
the direction of pi in Figure 3, the final configuration of the landmark points (say, after one unit 
of time) will differ greatly from the one caused by the original value of pi. Also, if you imagine q^ 
moving along two nearby parallel straight lines, the differential effect on the cloud of other points 
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accumulates so that the final configurations will differ everywhere; thus, even though the initial 
landmark configurations are close, the final configurations will be far away. In general, the last 
negative term in the curvature expresses the same effect: the global drag effect of each point results 
in a kind of turbulent mixing of all the other points (think of a kitchen mixer the motion of whose 
blades mixes the whole bowl). 

Proposition 14 simplifies in the case of £ = £^(]R^) (two landmarks only). We shall write: 

af :=(ai,ui2), := ai - af W12, := (/3i,«i2), 13^ := pi - pi U12. (57) 

Proposition 15. In the case of C = £^(M^), when a,(3 £ {T*£)i the numerator and denominator 
of sectional curvature are given by, respectively: 

R{a\pKpKai) = R, = -^70 \\PW " 4^"^ (58) 

4 70 + 712 

Mt'C Mt'C - {cx,P)t'C = 7o(||M'"^ - «f/3^f + kW®o.i - ai^Pif). (59) 

Proof It is the case that R4, = -|||iJ2||^(K~i)22. But (K~^)22 = (7o - 7i2)~Wo, whereas from 

Proposition 14 wc have 

||^^2|P = (70 -7i2)'((ai, VXi2)2||^^||2 + vO'll"i|l' - 2(ai, Vifi2)(/3,, vxi2)(ai,/?i)), 

where \/K^^ = ')i2U\2 by (53). Inserting expressions (57) into the above formula yields (58). 
Prom Proposition 11 we have that the denominator is given by 7o(||cki||^||/^i|P ~ (cHi A)^); again, 
inserting (57) into such formula yields (59). □ 

We will generalize the above results in the next section. 

5 Landmark geometry with two nonzero momenta 

The complexity of the formula for curvature reflects a real complexity in the geometry of the land- 
mark space. But there is one case in which the geometry of landmark space can be analyzed quite 
completely. This is when there are only two nonzero momenta along a geodesic. To put this in 
context, we first introduce a basic structural relation between landmark spaces. 

5.1 Submersions between landmark spaces 

Instead of labeling the landmarks as 1, 2, • • • ,N, one can use any finite index set A and label the 

landmarks as with a G A. And instead of calling the landmark space we can call it 
Now suppose we have a subset B <Z A. Then there is a natural projection tt : — > C'^ gotten by 
forgetting about the points with labels in A — B. In the metrics we have been discussing this is a 
submersion. In fact, the kernel of dn, the vertical subspace of T£-^, is the space of vectors v"' such 
that I)" = if a e B. Its perpendicular in T* is: 

(T*£^)b := {p G T*C^ \pa = for a G A- B} 

so the orthogonal complement of ker((i7r) in T£^ is the space of vectors p" where p is in (T*£-^)b. 
On this subspace, the norm is just 

b,b'el3 
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whether is taken to be a tangent vector to A or to B. In other words, the horizontal subspace 
for the submersion tt is the subbundle (T*£-^)g C TC-^ of tangent vectors where p has zero 
components in A — B and this has the same metric as the tangent space to JC^. In particular, 
from the general theory of submersions, we know that every geodesic in C'^ beginning at some 
point Tf{{q°'}) has a unique lift to a horizontal geodesic in starting at {g"}. The picture to have 
is that all the landmark spaces form a sort of inverse system of spaces whose inverse limit is the 
group of diffeomorpisms of . 

We don't want to pursue this is in general, but rather we will study the special case where the 
cardinality of B is two. We might as well, then, go back to our former t(;rniinology and consider the 
map TT : £^ -> gotten by mapping an iV-tuple {q^,q^, ■ ■ ■ , q^) to the pair {q^,q^). Moreover, we 
want to consider only the case in which the kernel K is rotationally invariant as in (52). A basic 
quantity in all that follows is the distance p := — between the two momentum bearing points. 



5.2 Two momentum geodesies 



Remarkably, we can describe, more or less explicitly, all the geodesies which arise as horizontal 
lifts from this map. These are the geodesies with nonzero momenta only at and g^. Moreover, 
the formula for sectional curvature for the 2-plane spanned by any two horizontal vectors can be 
analyzed. This analysis was started in the PhD thesis of the first author [19] and has been pursued 
further in [18]. 

The metric tensor oi C = C'^{R^) in coordinates is obtained by inverting the 2x2 matrix K; 



K 



70 l{p) 
l{p) 70 



(K-i)i2 = (K-i)2i = -(7o' - i{p?)-Hp) 



(60) 



so that the cometric and metric, for all covectors a, /3 e T*C and vectors v,w G TqC, are simply: 



g-\a,p) = 7o((ai, A) + {a2,p2)) +l{p){{ai,p2) + (a2,/3i)), 



(61) 



g{v,w) 



l{pf 



The geometry of the two-point space is best understood by changing variables for the landmark 
coordinates {q^,q^) and the momentum {p^jp"^) to their means and semi-differences, that is: 



q 



1+9 X 

, Sq: 

so that: q^ =q + 5q, q 
Then the cometric (61) becomes: 

g-\{a,5a), {^,5^)) = 2(70 +7(p)) <S,^) + 2(70 - 7(p)) (<5a, 



^ q - q 
2 

-2 = g - 



p :-- 



p^ +p^ 



p^ = p + Sp, P^ = p — Sp. 



5p :- 



P -P 
2 



(62) 



With these coordinates, the two-point landmark space becomes a product V x Vs in which all fibres 
V X {6qo} are flat Euclidean spaces though with variable scales, all fibres {^q} x Vs are conformally 
flat metrics sitting on the manifold M.^ — {0} and the tangent spaces of the two factors are orthogonal. 

Proposition 16. In terms of means and semi-differences, the geodesic equations for C'^{EP) are: 



'q={lo+ 7(P)) P, 

Sq = (70 - lip)) 5p, 



p = 0, 

sp--2'-^mf 



\dpf)dq. 



(63) 



24 



The above result is proven by direct computation. We can solve these equations in four steps. 

1. First the linear momentum p is a constant, so "center of mass" q moves in a straight line 
parallel to this constant: 

q{t) = m + ( ^ (70 + 7(/>(r)))rfr) p. (64) 

2. Secondly, if we treat vectors Sq and Sp as 1-forms in M^, equations (63) also show that: 

{Sq A Sp)' = Sq A 6p + Sq A Sp = [(scalar) • Sp] A Sp + Sq A [(scalar) • Sq] = 0, 

so the angular momentum, 2-form Sq ASp e /\^ M.^ is constant; we write this as w A where w 
is the nonnegative real magnitude of the angular momentum and (e^,e^) is an orthonormal pair. 
Then it follows that: 

5q{t) = y{t)[cos {6{t))e^ + sin {e{t))e^] , for some function e{t). 

3. Thirdly, we can express 6{t) as an integral: 

Sq = lp[cos{e)e^ + sin(6')e^] + - sin(6i)e^ + cos(6')e^] , so 
Sq ASq = - \p'^Oe^ A e^, as well as (from (63)): 
SqASq= (70 - 'y{p))SpASq = -^(70 - 7(p))e^ A e^; 
combining the second and third lines, we find: 

0{t) = e{0) + 4ujJ^ ^° ~ ^2^^^*^ (65) 

note that 9{t) is a m,onotone increasing fimction if cj 7^ 0, otherwise it is a constant. 

4. The last step is to solve for p{t) . This can be done using conservation of energy. Equations (63) 
are in fact the cogeodesic equations for the Hamiltonian 7i{p, q) of section 4.1, which we may rewrite 
in terms of means and semi-differences as 

n = {io + i{p))\\vf + bo-i{p))Wf 

by (62); hence this function of p and \\Sp\^ is a constant (p is also a constant). Then we calculate: 

{p^y =A{Sq,Sqy = m,Sq) = S{^n-^{p)){SpM) =^ p = A^-^^^^{SpM)■ 



But: 

{Sp, Sqf + = {Sp, Sqf + \\Sp A Sqf = \\5pf • = ^ - (7o + 7(p) )Hp"^ 

4 V 7o - lip) 

P 

This means that the function p{t) is the solution of: 

t= :r7=r^, where: F{x) := H x^^o - 7{x)) - Wp^x' {j'o - li^f) - {lo - lix))' . 

(66) 
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Figure 4: Converging and diverging trajectories for two landmarks in two dimensions. In these ex- 
amples ^(x) = exp(-ix2), (gi(0), g2(oj) = ((1, 0), (-1, 0)), (pi(0),p2(0)) = ((-10, 8.6), (10, -8.6)) 
for the graph on the left, (331(6), 332(0)) = ((-10,9), (10, -9)) for the graph on the right. The thick 
and thin lines are, respectively, the trajectories of the first and second landmarks. 

Summary. If wc fix constants "H, p, w, p(0), 6'(0), <?''(0) (for all a), we can first integrate (66) 
to get p(t) (the separation of and g^), then integrate (65) to find their relative angle ^(f), then 
integrate (64) to get their center of mass g(i). This gives the trajectories of and <f . The remaining 
points are dragged along as solutions of: ^q°'{t) = j{\\q°'{t) - q^{t)\\)pi{t) +j{\\q'^{t) — q'^{t)\\)p2{t). 

As worked out in [18, 19], one can classify the global behavior of these geodesies into two types. 
One is the scattering type in which q^,q^ diverge from each other as time goes to cither ±00. This 
occurs if the linear or angular momentum is large enough compared to the energy. In the other case 
where the energy is large enough compared to both momenta, they come together asymptotically 
at cither t = +00 or —00, diverging at the other limit. In both cases, they may spiral around each 
other an arbitrarily large number of times (see Figure 4). 

5.3 Decomposing curvature 

Next we consider £^(M^): we want to compute the sectional curvature R{a^, /?", a") for cotangent 

vectors that are nonzero at only (g^,?^). Also, we will use the notation u := ||^i~^2|| for the unit 

vector from to q^ as well as p = — for their distance. Similarly to (55), we will also want 
to decompose any vector in r/ e into its parts tangent to u and perpendicular to u: 

77" := (77, u), and rj^ := rj — ry" u. 

Once again note that t?" is a scalar whereas -q-^ is a vector. Following the notation used to describe 
geodesies above, for any a G {T*C)\^2 '■= {rj G T*C | rya = for o > 2}, we write a = \{a\ + a-2) 
and 5a = \{ct\ — 0:2) • 
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Proposition 17. In >C^(R^) for any pair a,(5 & {T*C)i^2, the terms R\,R2 and Rs in the numer- 
ator of sectional curvatMre can he written as 

i?i = 4(70 - 7(p)) V(/5) {Sj/3i - 6pKi,6J/32 - 6/3^2) 

+ 4(70 - 7(p))' ^ (Sa^ A - Sp^ ai,Sa^ ®p2- (E) as), 

R2 = -4(70 - lip)) i{pf <^all/3i - 5a^2 - ^/3"a2), 

i?3 = IfI_My(p)2[((«^,^2) + (/3i,a2))' -4(ai,a2)(/3i,/32)]. 

Once again we have used: (ui ®wi,V2® W2) ■= {vi,V2) {wi,W2) ■ We need the following result. 
Lemma 18. For any a € (T*£)i_2, the discrete strain Si2{a) is given by: 

5i2(a)=2(7o-7(p))^a- (67) 
For any pair a,j3 & (T*C)i^2 it is the case that Fa{a, /3) = for a > 2, whereas 

F,{a,l3) = -F2{a,0) = ^{{a,,P2) + {Pi,a2))u. (68) 

Also, 

D\a, (3) = 2(70 - lip)) iip) 5 J p2. D^{a, /3) = 2(70 - 7(p)) l\p) ^a" /^i- (69) 

Remark. We are not interested in D°'{a,[3) for a > 2 since the terms in formula (42) where they 
appear are zero (because Fa{a, /3) = for a > 2). 

Proof of Lemma 18. The formula for the discrete strain results from: 

Si2{a) = (a«)i - {a^f = Ebi^^^ ' K^^hb = 70^1 + 7(^)^2 - 7(p)ai - 70^2 = 2(70 - l{p))Sa. 
The values for F follow immediately from formula (38) and VK^'^ = Y{p) u. Note that: 

Ci2(a) = C2i{a) = (5i2(a), VXi2) = (2(7^ - j{p))Sa,i'ip)u) = 2{jo - l{p))l' (p) SJ , 
and similarly we have that Ci2(/3) = C2i(/3) = 2(70 - 7(p))7'(p) SP^. So 

D\a,p) = Ci2ia)p2 = 2(70 - 7ip)h'{p)Sa^2, 

D\a,p) = C2i(a)A = 2(70 - 7(p))7'(p)^allA- □ 

Proof of Proposition 1 7. The Ri expression follows by substituting the expressions in (67) into for- 
mula (56), noting that the only non-zero terms in the latter are for (a, 6) = (1, 2) and (a, 6) = (2,1). 
By Theorem 9 and the fact that F2 = — Fi from Lemma 18, R2 is given by 

R2 = {D\a, a) - D\a, a),F^{P, /?)) + {D\p, /3) - D\(}, /J), Fi(a, a)) 

- {D\a, p) - D\a, p) + D\p, a) - D\p, a), Fi(a, p)). (70) 

Again by Lemma 18 we have that D^irj^C,) — D'^{r]X) = —4(70 — 7(p))7'(p) S( for any pair 

r?, C e (r*£)i,2, while Fi{r],C) = 57'(p)((»?i> C2) + {'n2,Ci))u- Applying this to all the terms we get 
the expression for R2 in the statement of the proposition. 
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As far as Rs is concerned, by Theorem 9: 

R3=7o[{Fi{a,l3),Fiia,l3))-{F,{a,a),F,{0,0))]+j{p)[{F,{a,l3),F2{a,l3))- 

+ 7(p) [(F2(a, /3),i^i(a, /?)) - (i^2(a, /?))] +7o [(J^2(a, /3),F2(a, - (F2(a, a),F2(/3, 

= 2(70 - 7(p)) [(i^i (a, i^i («. 0)) - {Fi {a, a),F, (/?, /?))] 



7o-7(p) 
2 



(Y [((ai, ^2) + Q2»' - 4(ai, a2){f3i, P^)] , 



where we have used the fact that i^2 = — -F^i, by equation (68). This completes the proof. □ 

The expressions provided by Proposition 17 become much clearer if we go over to means and 
semi-differences, i.e. if wc use the substitutions: 

ai = a + 5a, a2=a-Sa, (3i=]3 + 5/3, /Ja = ^ - 5(3. (71) 

Corollary 19. For any a,/3 e (T*£)i,2, with C = £^(M^), it is the case that: 

iJi = 4(70 - 7(/9))%"(/5) {\\5l3\\a - 5a\\pf - \\5/3Ua - Sahpf) 

+ 4(70 - -f{p)f'^^{\\5l3^(g)a- 5a^ ^pf - j|(5/3^ (g> 5a - 5a^ 5pf), 
i?2 = -4(70 - 7(p))7'(/o)' - 5JW " - 5J5(3\\^), 

i?3 = {jo-l{p))l'{p? {2\\5l3(g)a-5a(g}W - (E)a -a<g)P\\^ - \\5j3 <g) 5a - 5a<g) 5p\\^) . 
Proof. By insertion of formulae (71) it is easily seen that 

{6a^i - 5^11 ai, foil /32 - 5^lla2) = \\5(3^\a-5a^f - \\5pha - 5^5/3^, 
{5a^ OA- 5p-^ (g) ai,5a-^ O /32 - 5p-^ O 02) = \\5l3-^ (g) 5 - fo-^ O pf- \\5l3-^ (g>5a- 5a-^ O 5/3 f, 

so the new expressions for Ri and R2 follow immediately. Also: 

[(ai , A) + , a2)] ' - 4(ai , as) (A , A) 

= [2{{a, p) - {5a, 5(3))] ' - 4((5, a) - {5a, 5a)) {{p, p) - {5p, 5p)) 
= -2[2{{a,a){p,p) - {a,p)^)] - 2[2{{5a,5a){5p,5p) - {5a, 5p)^)] 

+ 4[{a,a){5p,5p) + {p,p){5a,5a) - 2{a,p){5a,5p)] 
= -2\\P O a - 5 O P\f - 2\\5P ^5a-5a^ 5Pf + 4\\5P 5 - 5a (g) P\f. □ 

The fourth term R4 is the only one which involves the other points q"',a > 2. But one has an 
inequality for this term involving the same expressions in a and /3: 

Proposition 20. Any pair a,p £ {T*£^)i^2 o-^e constant 1-forms on £^ which are pull-backs via 
the submersion — >■ of constant 1-forms on C? . We can therefore consider the curvature term 
RAiC^) = -fllM, IP on and the corresponding term Ri^C^) = -l\\[oiKP%2f on C? . 
Then we have the inequality: 



RAiC"") < R4C^) = -&i{pf 



(70 7(p)) ||g^||^_^^||^||2^ / \\sp\\5a-5a\ 
70 + 7(P) 
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„1 „2 



^02 



ai a2 



Figure 5: Typical covectors a = (ai,a2) in spaces JUt*/;, (J-^T*/;, and T*/;, for £ = /;2(k2)_ 

Proof. Firstly, note that [a'j^'J^w breaks into perpendicular parts: a vertical part in the kernel of 

diT and a horizontal part which is simply the horizontal lift of [a\ (i'^]c'^ ■ This explains the inequality 
assertion in Proposition 20. To calculate i?4(£^), we use the last expression in (44), i.e. 



70 + 7lPj L 

- 27(p)(fcll/3i - ^^llai,(5all/32 - ^/3"a2)| 



where we have used (60) and (69). The final result follows after inserting (71) into the above 
expression and performing some algebra. □ 

Note that all terms in Corollary 19 and Proposition 20 are very similar. In fact, they are all 
"components" of the norm Ha A of the 2-form whose sectional curvature is being computed. 
First note that we can decompose T* into the direct sum of three pieces, namely: 



|(au, — a?i) I a G M}, 

{(p, -p)|peM^,p±M}, 
{(p,p)|peM^}, 



dim((5llT*£2) = 1^ 
dim(5-Lr*£2) =£)-!, 
dim(r*£2) 



where as usual u := n^i^^ay (see Figure 5). Note that these throe subspaces are orthogonal with 
respect to the cometric by virtue of (61). An arbitrary covector a = (ai, a2) G T*C^ can be uniquely 
decomposed into the summation a = a(i) + a(2) + Q^(3)) with: 

«(!) := {SJu, -SJu) G S\^T*C^, a(2) := (fo-^, -Sa-^) G 5^T*C'^, a^) := (a, a) e T* C"^ . (72) 

So it is the case that: (i) a € (5llr*£2 ^ J-Lq, = and 5 = 0; (ii) a G J-Lr*/:^ ^ = and a = 0; 
(iii) a e T*£2 ^q,|| = q and ^a-^ = 0. 

Consequently the space of 2-forms /\^ T*£^ decomposes into the direct sum of five pieces: 



/\r*£2 = 0y,, with: 

i=l 

^2 := ^^T*£2 A T*£ 

2 

Vi := /\{S^T*£^), 



Vi := (5"T*£2 a r*£2, 
:= 5llT*£2 A ^-LT*£2, 

^5 := Ar-c')- 
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(Since (5llT*£2 

is one-dimensional it creates no 2-forms.) Once again, note that the spaces Vi, . . . , V5 
are pairwise orthogonal with respect to the inner product 

(q A ^, ^ A ri)^2 := {a, Or^r^ ('^' ^)t*£2 " ^)t'C^ Ot-c^ ' /5> C> e T*C'' (73) 

by the orthogonality of 6^^T*£., 5^T*C, and T*£. Any 2-form a A /3 then decomposes into the sum 
of its five projections onto these subspaces and its norm squared is the sum of the norm squared of 
these components. Let us first give the five pieces of its norm names: 

Ti := ||(5^"u(^a- (5a"u(g)/3||^ 
T2 := \\6I3^ ®a-5a^ T3 := ||(5/3llu fo-^ - foUu (g) J/J-^f, 

Ti := \\5^^ ®5a^ -5a^ ®5p^\\^, T5 := ||^ a - a (g) ^||^ 

In the above definitions || j| indicates the Euclidean norm. We have to be careful here: we have 
been using Eiiclidcan norms in in all our formulas above and now we arc dealing with norms in 
T* C^-. these essentially differ only by a factor, by (62). More precisely, the following result holds: 

Proposition 21. The denominator of the sectional curvature (16) /or £^(R^) can he written as: 

\\a A = 4(72 - 7(p)')(Ti + T2) + 2(70 - 7(P))'(2T3 + T4) + 2(70 + ^{p))^n. (74) 

Proof. We may apply decomposition (72) to both a = X]f=i (^{i) ^^'^ P = /^(j)) and write 

aAP =(«(!) A /3(3) - A a(3)) + (a(2) A /3(3) - ^(2) A a(3)) 
+ ("(1) A /3(2) - A q;(2)) + a(2) A /3(2) + a(3) A ^(3), 

where the five summands on the right-hand side belong to Vi, . . . , V5 respectively. We have 

A /3(3) - A a(3)||^2y.£2 = 
= ||a(i) A /3(3)ll^2 j,.£2 + ||/3(i) A Q;(3)||^2 j,.£2 - 2(a(i) A /3(3),^(i) A a(3))A='T*£2 

= 4(7o' - lipf) [{Sa^mf + {SP^^rm' -26JSl3Ha,P)] = 4(70^ - ^ipf) T,, 

where we have used (73) and (62) in step (*). The square norm of the remaining four terms is 
computed similarly. Orthogonality of Fi, . . . , V5 finally yields (74). □ 

To express the formulas for the numerator of sectional curvature succinctly, let us also introduce 
abbreviations for the coefiicients involving 7: 

fci(p):=(7o-7(p))V(p), Hp) 

k3{p)--={lo-lip)h'{pf, h{p) 

Note that ki, fc2, ^3 and k4 are all homogeneous of degree 3 in 7 and degree —2 in the distance p 
or dp on £^ . Moreover ^2 is negative, ^3 and ki are positive, while fci may be positive or negative. 
For all 7 of interest, 7' is everywhere negative, starting at decreasing to a minimum at some po, 
then increasing back to at 00. Then ki is negative for p < po and positive for p> po. 
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(70 - 7(P)) 
(70 -7(P))' 



27 [p) 



(75) 



70 + 7(P) 



I'ipf- 



Figure 6: The coefRcients of Ti (top left), T2 (top right), T3 (bottom left) and T4 (bottom right) for 
the Bessel kernels 7 (shown with thin lines) and the Gaussian kernel (shown with the thick line). 
The kernels are scaled to normalize 7(0) and 7"(0). 

The following equalities are proven by direct computation: 

\\Sp^a-Sa^W = Ti+T2, 
\\6l3-^ (g)Sa- 5a^ (g) Sl3\\^ + T4, 
\\SP (S)6a-Sa® Sf3\\^ = 2T3 + T4. 

Inserting notation (75) and the above equalities into Propositions 17 and 20 immediately yields: 
Proposition 22. We can write the terms in the numerator of sectional curvature for C^{R^) as: 

i?l =4fci(Ti-T3) + 4fc2(T2-r3-T4), i?2 = -4fc3(Ti-T3), 

R3=k3{2{Ti+T2)-2T3-Ti-Ts), = -Qik^T^ + k^T^), ^ ' 

hence R = R{aK , a«) = J2Ui ^ 

may he expressed as: 

R = 2(2fci - fcs - 3fc4)Ti + 2(2fc2 + fc3)T2 + 4( - fci - fca - fc3)T3 + ( - 4^2 - ^3)74 - fcsTs. (77) 

By virtue of Proposition 20 the above proposition still holds in the case of £^{M.^) as long 
as a, /3 e (T*£)i_2 and the equality signs for R4 in (76) and R in (77) are substituted by "<". The 
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GEODESIC 1 

GEODESIC 2 

-0 2 I ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' ' 

0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 2 4 6 8 10 12 

P q' 

Figure 7: Left: sectional curvature JC for £^(M'^) (from Proposition 23), as a function oi p = Iq^ — q^l; 
here 7(0;) = exp(— ^x^). Right: two trajectories in £^(M^) shown in the {q^,q^) plane (under the 
assumption that < q^). Both geodesies originate at [q^^q^] — (0,4), and lie in the region 
where /C > (above the upper dotted line, that indicates the zero of K. at \q^ — q^| ~ 1.53); they 
have different initial momenta {pi,P2) = (1, 1) and {pi,P2) — (1,0.4), and exhibit conjugate points. 

coefficients in (77) may have all sorts of signs for peculiar kernels. However, the kernels 7 of interest 
are the Bessel kernels (3) and the Gaussian kernel, which is their asymptotic limit as their order 
goes to infinity. The coefficients for these kernels are shown in Figure 6. We see that the coefficients 
of T2 and T3 are negative while those of T-i are positive. Henceforth, we assume we have a kernel 
for which this is true. 

5.4 Sectional curvature of £^(]R^) 

Finally, we will now explore the important example of two landmarks on the real line. In this partic- 
ular case the manifold is two dimensional, so sectional curvature /C will turn out to be independent 
of cotangent vectors a and /3. In fact, given the translation invariance of the metric tensor, it will 
only depend on the distance p = \q^ ^ q^\ between the two landmarks. 

The spaces S'^^T* C'^{M}) and T* C'^{M}) are one-dimensional while 5^T*C:^{M}) = {0}. Thus 

2 

f\T*C^{M}) = 5llT*£2(Mi) A T*C^{M}) 

and the only non-zero term in (77) is Ti. Therefore combining formulas (74) and (77) we get: 
Proposition 23. The sectional curvature of C'^iJS}) is given by 

^ 2ki - fcs - 3fc4 ^ 7o - lip) II _ 27o~7(p) , .2 
2(7o' - lip?) 70 + lip) (70 + lip))' ■ 
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■ - GEODESIC 1 
— GEODESIC 2 



Figure 8: Existence of conjugate points in £^(K^), with j{x) = exp(— Both geodesies 
originate at landmark set {q^,q^) = ((— 2, — 4), (2, 4)) ; the first one (dashed) has initial momen- 
tum {pi,P2) = ((0, fO), (0, 10)) e r*£^ while the second one (continuous) has initial momentum 
{pi,P2) = ((6, 10), (-6, 10)) e 5llT*£2 e T*£^. The geodesic trajectories exhibit conjugate points. 



The above function /C is shown on the left-hand side of Figure 7 as a function of p, for the 
Gaussian kernel. The coefficient of the term Ti in (77) is negative for p small and positive for p 
large. The "cause" of the positive curvature has been analyzed in [19]. Roughly speaking, suppose 
two points both want to move a fixed distance to the right. Then if they are far enough away, they 
can just move more or less independently (we shall refer to this as Geodesic 1). Or (i) the one in 
back can speed up while the one in front slows down, then (ii) when the pair are close, they move in 
tandem using less energy because they are close and finally (iii) the back one slows down, the front 
one speeds up when they near their destinations (Geodesic 2). This gives explicit conjugate points 
and is illustrated on the right-hand side of figure Figure 7 (where Geodesies 1 and 2 are represented, 
respectively, by the dashed and thick curves). 

5.5 Sources of positive curvature; obstacle avoidance 

There is another source of positive curvature in in higher dimensions. It is clear from equation (77) 
and Figure 6 that any positive curvature must come from the term with Ti or the term with T4. As 
the five terms are orthogonal, we can make all of them but one zero. 

For example, if we choose a = (da^^u, —Sa^^u) <E (5"T*£ and /? = (/3,/3) € T C, then it is the case 
that Ti = (i5q;II)^||/3|P and it is the only non-zero term. Then, if p is sufficiently large, the sectional 
curvature for this 2-plane is positive as discussed in the last section. Figure 8 illustrates an instance 
of the existence of conjugate points for two geodesies in £^(R^); the momenta (pi,P2) of each of the 
two trajectories belong at all times to d^^T*jC'^ T*£^. 

The other possibility is that T4 is the non-zero term, which happens when a = {6a-^, —Sa-^) G 
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5^T*C and /3 = (5p^. -Sd^) e S^T*C. Wc have T4 = 2{\\Sa-^f\\dl3-^\\'^ - {Sa^ ,5^-^)'^), and for it 
to be nonzero it is required that D > 3 because T4 is the norm of a 2-form in /\^ {5^T*C), which 
has dimension {D —1){D — 2)/2. The positive curvature of this section is readily seen by considering 
the geodesies which these vectors generate. The simplest example is the following: 

Proposition 24. The circular periodic orbit of radius r: 

(t) = (r cos t, r sin t) , {t) = -q^ {t) , (78) 

t gM., is a geodesic in if and only ifr is the solution of the equation 70 — 7(.t) + 2x'y'{x) = 0. 

Proof. For orbit (78) it is the case that q = 0, Sq = 0, p = 2r, p = 0, Sp{t) = (70 — 7(p)) ^q^{t)l 
therefore equations (63) are satisfied if and only if 70 — 7(r) + 2r'y'{r) = 0. □ 

(The above result was also proved by Frangois-Xavier Vialard of Imperial College, London.) 
Orbit (78) has the property that at time tt, q^ and interchange their positions: it is a geodesic 
from the set of landmark points ((r, 0), (-r, 0)) G C^iR"^) to the set ((-r, 0), (r, 0)) G C^iR"^). But 
if these points live in K^, they can move around each other in any plane containing the points. Thus 
we have a circle of geodesies in £^(R^): 

q^ (t) = (r cos t, r cos 6 sin t, r sin 6 sin t) , q^ [t] = —q^ (t) 

all connecting ((r, 0, 0), (-r, 0, 0)) to ((-r, 0, 0), (r, 0, 0)) , for any 6 G [0, 27r). This is exactly like all 
the lines of fixed longitude connecting the north and south pole on the 2-sphere and means that one 
set of landmark points is a conjugate point of the other in £^(M^). This is the simplest example of 
how geodesies between landmark points must avoid collisions and so make a choice between different 
possible detours, leading to conjugate points and thus positive curvature. 

6 Conclusions 

We have computed a formula for sectional curvature for C'^(R^), the Riemannian manifold of N 
landmark points in D dimensions. To do so we have developed a formula to compiite sectional 
curvature of a Riemannian manifold in terms of the cometric, its partial derivatives, and the; metric 
(but not its derivatives). Finally, we have fully examined the case of geodesies in which only two 
points have non-zero momenta, and found that there are essentially two sources of positive curvature; 
one only occurs when D > 3. Future work may include: exploring the shape of the coefficients in (77) 
for different kernels; finding new sources of positive curvature when momenta are non-zero at more 
than two points; analyzing what happens asymptotically when the points are very close or very far 
from each other; further relating the dynamics of landmarks to the geometry of the space. 
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